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INTRODUCTION  t£{- 99 ~  ^ 

The  flow  of  vibrational  power  through  structures  is  an  important  characteristic  in 
understanding  the  behavior  and  the  dynamic  response  of  an  elastic  structure.  In  the  case  of  thick 
or  curved  structural  components,  the  vibradonal  power  can  be  transmitted  by  different  wave  types 
which  have  different  propagation  characteristics  through  the  structure  and  which  can  propagate  at 
different  speeds.  A  wave  propagating  through  a  structure  can  be  scattered  into  other  wave  types 
when  the  wave  encounters  discontinuities  in  the  structure,  such  as  changes  in  cross  section,  stiffeners 
or  bulkheads,  etc.  If  the  structure  is  fluid  loaded,  some  of  these  waves  propagate  at  supersonic  speeds 
and  can  therefore  very  efficiently  loose  some  of  their  energy  to  acoustic  radiation.  The  efficiency 
of  noise  radiation  from  a  structure  is  very  much  dependent  on  the  presence  of  wavenumber 
components  which  are  supersonic  relative  to  the  acoustic  fluid  that  surrounds  the  structure.  The 
supersonic  wavenumber  components  can  be  associated  with  either  waves  that  propagate  at 
supersonic  speeds  or  with  supersonic  wavenumber  components  created  due  to  the  presence  of 
discontinuities.  Consequently,  for  noise  control  purposes,  it  would  be  required  to  differentiate 
between  the  relative  contribution  and  strength  of  the  structural  intensity  that  can  be  associated  with 
different  wave  types  or  wavenumbers. 

Additionally,  dissipation  mechanisms  introduced  to  dissipate  the  propagating  vibrational 
power  can  be  wave  type  dependent.  Dissipation  efficiency  is  different  according  to  the  type  of  wave 
and  certain  type  of  waves  can  act  as  short  circuiting  paths  for  the  energy  transmission.  Therefore, 
an  understanding  of  the  power  propagated  by  different  wave  components  and  how  the  power 
propagated  by  one  wave  type  can  be  scattered  into  other  wave  types,  is  an  essential  element  in  the 
design  for  vibration  control.  A  formulation  is  required  to  relate  the  scattering  of  vibrational  power 
from  one  wave  type  to  another  based  on  the  characteristics  of  the  discontinuity.  The  control  of 
scattering  of  vibrational  power  into  those  wave  components  which  are  either  efficient  radiators  or 
can  act  as  short  circuits  around  energy  dissipation  mechanisms  can  reduce  the  sound  radiated  by  a 
structural  component 


The  objectives  of  this  work,  as  defined  at  the  beginning  of  the  project  were  to  develop 
analytical  and  experimental  techniques  for  the  analysis  of  power  flow  through  thick  connected  plate 
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structures  and  to  identify  the  components  of  vibrational  power  flow  that  can  propagate  across  a 
discontinuity  in  a  thick  plate  structure.  The  emphasis  was  to  be  placed  on  the  interaction  and 
scattering  between  the  different  wave  components  in  the  presence  of  structural  discontinuities.  An 
approach  was  to  be  developed  that  can  directly  provide  the  structural  intensity  components 
associated  with  the  different  wave  types.  The  influence  of  fluid  loading  was  to  be  considered  based 
on  the  progress  of  the  research.  Specifically  the  objectives  were  as  follows: 

1.  Using  the  formulation  for  analytically  describing  the  interaction  between  in-plane  and 
out-of-plane  waves  developed  in  the  previous  work,  formulate  the  basis  for  determining  the 
influence  of  discontinuity  characteristics  on  the  scatter  of  vibrational  power  across  the 
discontinuity. 

2.  Experimentally  measure  the  components  of  the  structural  intensity  associated  with  the  different 
wave  components. 

3.  Develop  the  formulation  for  integrating  the  influence  of  fluid  loading,  and  the  additional 
coupling  that  exists  through  the  fluid  for  coupled  structures  for  finite  coupled  plates. 

At  the  beginning  of  this  research  program,  wo  techniques  were  available,  in  the  open 
literature  that  can  be  used  to  evaluate  the  flow  of  vibrational  power  through  thick  connected 
structures,  especially  if  tire  structure  is  finite  and  mode  dependent  results  are  required.  Some  work 
had  been  done  previously  for  two  plates  joined  in  an  L-shaped  configuration  [  1]  which  did  take  into 
account  both  the  in-plane  and  the  out-of-plane  waves.  However  in  that  analysis  [1],  a  thin  plate 
formulation  was  used  for  the  representation  of  the  in-plane  waves  and  the  coupling  between  the 
in-plane  and  out-of-plane  waves  only  existed  because  of  the  configuration  of  the  structure  being 
analyzed. 

Furthermore,  the  concept  of  structural  intensity  or  power  flow  had  not  been  fully  explored 
as  a  method  by  which  one  can  describe  the  behavior  of  elastic  structures.  The  type  of  information 
that  can  be  obtained  from  structural  intensity  analysis  includes: 

1.  Location  of  vibration  sources  and  sinks 

2.  Critical  power  flow  paths 

3.  Locations  on  the  structure  with  a  high  level  of  energy  loss  associated  with  either  acoustic 
radiation  or  mechanical  dissipation. 

4.  Identification  of  the  dominant  wave  components  (with  a  component  decomposition  of  the 
structural  intensity)  that  are  involved  in  the  mechanism  by  which  the  vibrational  power  is 
transmitted. 
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had  not  been  considered. 


With  regards  to  experimental  analysis  of  the  power  flow,  past  work  generally  has  been 
limited  to  thin  plates  or  beams  of  infinite  extent.  Some  work  has  been  done  by  Williams  [2]  that 
considers  thin  elastic  shells  and  which  takes  into  account  propagation  by  both  in-plane  and 
out-of-plane  waves.  The  work  by  Pavic  [3]  on  thick  structures  concentrates  on  developing  a 
measurement  techniques  using  strain  gauges  and  velocity  transducers  to  measure  the  in-plane 
velocities.  While  this  paper  contains  a  very  good  discussion  on  the  errors  associated  with  the 
measurements,  the  approach  is  not  easy  to  implement  because  of  the  required  installation  of  the 
strain  gauges  and  the  way  the  in-plane  motions  are  being  detected.  Berthelot  and  Jarzynski  [4] 
designed  a  laser  probe  for  the  measurement  of  in-plane  structural  motion.  Using  two  interfering 
laser  beams  the  out-of-plane  motion  of  the  structure  is  filtered  out  of  the  measurement  Both 
frequency  domain  results  and  time  domain  results  have  been  obtained.  By  scanning  the  laser  probe, 
the  in-plane  structural  intensity  can  be  determined.  While  this  is  one  of  the  most  useful  techniques 
for  measuring  in-plane  motion  or  structural  intensity,  different  probes  or  processing  would  be 
required  to  measure  simultaneously  the  out-of-plane  structural  intensity.  It  would  be  desirable  to 
have  a  single  probe  that  can  simultaneously  measure  both  the  in-plane  and  the  out-of-plane  wave 
motion  and  thus  structural  intensity.  None  of  the  available  experimental  technic 
between  structural  intensity  components  propagated  by  different  wave  types. 

The  work  done  can  be  classified  under  three  categories: 

1.  Analytical 

2.  Simulated  and  direct  experiments 

3.  Signal  Processing. 

ANALYTICAL  WORK 
Thick  Plate  Analysis 

The  analytic  approach  used  in  the  research  under  this  grant  is  based  on 
Flow  (MPF)  method  [5].  In  this  approach  a  complex  structural  element  is  divided  into  subsystems 
where  each  subsystem  represents  either  a  physical  subcomponent  of  the  global  structure  or  a  wave 
component  within  a  physical  substructure.  Mechanical  mobility  functr  ns  are  then  derived  for  each 
of  the  subsystems  without  reference  to  the  other  subsystems  or  the  type  of  loading.  Assumptions 
can  be  introduced  on  the  junction  loading  to  simplify  the  problem  by  reducing  the  number  of 
mobility  functions  included  in  the  analysis,  the  response  of  the  global  structure  is  then  obtained  by 
developing  a  mobility  matrix  equation  based  on  continuity  of  forces  and  displacements  (velocity) 
at  the  junctions  between  the  subsystems.  This  approach  has  the  advantage  that  by  dividing  the  global 
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structure  into  subcomponents,  analytical  solutions  can  be  obtained  for  relative  complex  structural 
systems  which  would  be  difficult  to  analytically  model  otherwise.  Additionally,  since  the  mobility 
function  required  in  the  matrix  equations  are  generated  for  the  subsystems  independent  of  the 
influence  of  the  external  loading  or  the  other  systems,  the  MPF  approach  is  very  efficient  when 
performing  parametric  analysis  on  the  response  of  the  global  structure,  for  variation  of  one  of  the 
parameters  of  one  of  the  subsystems.  Only  the  mobility  functions  associated  with  that  subsystem 
needs  to  be  re-evaiuated. 

The  MPF  technique  is  particularly  applicable  in  the  middle  frequency  range,  where  the 
middle  frequency  range  is  defined  as  being  higher  than  the  first  one  or  two  modes  of  the  global 
structure  which  can  be  very  effectively  determined  using  finite  element  (FE)  analysis,  but  below  the 
region  where  the  modes  become  strongly  coupled  and  other  techniques  such  as  Statistical  Energy 
Analysis  (SEA)  would  be  more  appropriate.  Beyond  the  first  couple  of  modes  of  the  global 
structure,  the  computational  requirements  of  the  FE  approach  may  be  prohibitive  while  the  results 
generated  by  the  SEA  method  would  be  too  generic  possibly  leaving  out  some  important 
characteristic  of  the  response  of  the  global  structure. 

Using  the  MPF  method,  a  solution  has  been  found  for  the  determination  of  the  structural 
power  flow  between  two  coupled  finite  plates  of  different  thickness.  In  the  solution  coupling 
between  the  in-plane  vibrational  motion  or  power  flow  and  the  out-of-plane  loading  has  been 
achieved.  Initially  the  solution  was  based  on  the  Mindlin  [6]  and  Kane  and  Mindlin  [7]  solutions 
for  out-of-plane  and  in-plane  analysis,  respectively,  of  thick  plates,  and  by  defining  the 
out-of-plane  external  load  acting  on  one  of  the  surfaces  of  the  plate  as  being  made  up  of  two 
components,  acting  on  opposite  surfaces  of  the  plate.  These  two  forces  have  two  resultant  at  the 
mid-plane  of  the  plate  where  the  solutions  of  Mindlin  are  derived.  One  resultant  load  is  a  transverse 
force  and  the  second  resultant  load  is  a  ’’pinching”  [12]  force.  The  influence  of  this  decomposition 
is  that  the  out-of-plane  load  acting  on  one  of  the  surfaces  of  the  plate  imparts  vibrational  power  to 
both  the  in-plane  and  out-of-plane  waves.  Using  this  set  up,  a  solution  for  the  structural  power  flow 
is  obtained  for  a  number  of  different  thickness  combinations  of  two  coupled  finite  plates.  The 
thickness  variations  range  form  equal  thickness  to  a  thickness  of  the  receiving  plate,  which  is  the 
plate  with  no  external  excitation,  being  two  times  thicker  than  the  source  plate,  which  is  the  plate 
with  the  external  excitation.  The  thickness  of  the  source  plate  is  kept  constant  in  the  analysis. 

The  results  of  the  first  analysis  [8]  show  that  for  the  out-of-plane  flow  of  vibrational  power 
between  the  two  plates,  any  variations  in  thickness  will  reduce  the  magnitude  of  the  vibrational 
power  flow.  Maximum  vibrational  power  transfer  occurs  when  the  two  plate  components  are  of 
equal  thickness.  The  results  for  the  in-plane  power  flow  are  however  different.  In  this  case,  as  the 
receiver  plate  thickness  increases,  the  magnitude  of  the  vibrational  power  flow  increases,  reaching 
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a  maximum  for  a  thickness  ratio  of  between  2  and  4.  With  further  increase  in  the  thickness  ratio,  the 
flow  of  vibrational  power  between  the  two  plates  decreases.  The  increase  in  the  magnitude  of  the 
power  flow  is  more  significant  at  low  frequencies,  where  the  product  of  the  longitudinal 
wavenumber  and  thickness  (kh)  is  less  than  approximately  O.S.  Because  of  the  fact  that  the  flow  of 
vibrational  power  associated  with  the  in-plane  waves  increases  as  kh  increases  (for  kh  <  O.S),  and 
the  rate  of  decrease  of  the  power  flow  associated  with  the  in-plane  waves  for  kh  >  O.S  is  slower  than 
that  for  the  power  flow  components  associated  with  the  out-of-plane  waves,  more  of  the  available 
vibrational  power  goes  into  in-plane  waves  as  the  plate  thickness  changes.  Therefore  changes  in 
structural  thickness  to  control  the  flow  of  out-of-plane  waves  vibrational  power  may  result  in  more 
power  being  placed  in  the  in-plane  waves. 

From  initial  the  work  based  on  the  equations  of  Mindlin’s  [6]  and  Kane  and  Mindlin’s  [7], 
a  modified  version  has  been  developed  which  uses  the  approach  of  Mindlin  and  Medick  [9]  which 
includes  the  second  thickness  harmonic  for  the  in-plane  waves.  The  frequency  at  which  the  second 
harmonic  becomes  important  is  comparable  to  the  frequency  of  the  first  thickness  stretch  mode  for 
structures  with  a  material  which  has  a  value  of  poisson’s  ratio  of  order  0.3.  This  is  common  for  most 
materials  and  therefore  Mindlin  and  Medick ’s  approach  is  more  appropriate.  Based  on  this 
approach,  the  same  type  of  analysis  for  determining  the  coupling  between  the  out-of-plane  and  the 
in-plane  waves  have  been  used.  In  general  the  results  are  comparable  with  previous  results.  That 
is,  for  the  problem  of  the  two  coupled  plates  with  different  thicknesses,  the  ratio  of  the  in-plane  wave 
power  to  the  out-of-plane  wave  power  increases  as  the  thickness  of  the  plates  increases.  The  ratio 
of  the  transmitted  power  to  the  input  power  generally  varies  in  the  same  way  as  was  previously 
obtained. 

The  results  obtained  thus  far  are  for  a  discontinuity  between  two  coupled  finite  plates 
consisting  of  a  change  in  the  thickness  between  the  two  plates.  The  same  approach  can  be  extended 
to  other  types  of  discontinuities  in  a  rather  straight  forward  manner.  There  is  a  limitation  on  the 
maximum  frequency  or  thickness  that  can  be  considered  in  the  analysis,  if  the  Kane  and  Mindlin  [7] 
approach  is  used.  However,  this  limitation  is  removed  if  the  approach  by  Mindlin  and  Medick  [9] 
is  used. 


One  interesting  new  feature  which  was  observed  while  performing  the  analysis  using  the 
modified  approach,  but  which  is  not  a  result  of  this  new  approach  is  that  no  scattering  takes  place 
between  the  vibrational  power  of  one  wave  type  to  another  at  the  discontinuity.  Initially  this  was 
thought  to  be  a  problem  with  the  analysis  but  when  analyzing  the  governing  equations  of  motion 
and  the  results  more  closely  it  was  observed  that  scattering  would  only  take  place  when  a  load, 
external  to  the  plate  is  applied  at  the  discontinuity.  The  change  in  the  plate  thickness  discontinuity, 
because  of  the  fact  that  the  Mindlin  approaches  are  for  the  mid-plane  of  the  plates,  with  a  Legendre 
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polynomial  expansion  to  describe  the  thickness  dependency,  no  additional  loading  at  the  surface  of 
the  plates  is  induced  by  the  presence  of  the  discontinuity.  The  discontinuity  only  influences  the  flow 
of  vibration  power  by  one  wave  type  or  the  another.  Any  additional  loads  are  created  at  the 
mid-plane  of  the  two  plates  which  therefore  do  not  couple  the  in-plane  with  the  out-of-plane  waves. 
This  results  in  no  interaction  between  the  wave  types  at  the  discontinuity  and  hence  no  scattering 
of  vibrational  power.  For  scattering  to  take  place,  a  discontinuity  which  creates  an  external  loading 
on  the  surface  of  the  plates,  such  as  would  be  obtained  from  a  stiffener  or  a  bulkhead,  is  required. 
If  this  is  considered,  there  is  scattering  and  the  significance  of  the  scattering  increases  as  the  force 
or  moment  created  by  the  discontinuity  increases. 

The  force  or  moment  created  by  the  discontinuity  (stiffener  or  bulkhead)  is  a  function  of  the 
impedance  of  the  discontinuity  and  the  velocity  of  the  plate  along  the  line  of  attachment  of  the 
discontinuity.  A  complete  analysis  on  the  variation  of  the  strength  of  the  scattering  with  changes 
in  the  impedance  or  mobility  of  the  discontinuity  as  a  function  of  frequency  has  not  yet  been 
completed.  From  this  analysis  it  would  be  expected  that  the  optimum  impedance  or  mobility  for 
minimum  scattering  of  the  vibrational  power  from  the  out-of-plane  waves  to  the  in-plane  waves 
can  be  obtained. 

To  verify  the  approach  used  for  coupling  the  in-plane  and  out-of-plane  waves  when  a  load 
is  applied  to  the  surface  of  the  plate,  based  on  the  Legendre  polynomial  expansion,  the  results  from 
the  Mindlin  and  Mindlin  and  Medick  analysis  have  been  compared  to  results  obtained  from  a 
solution  of  the  three  dimensional  stress  equations  using  an  approach  similar  to  the  one  proposed  by 
Hutchinson  [10].  There  are  some  restrictions  when  using  the  3-D  stress  solution  without  the 
approximation  of  Mindlin  [6,9].  One  of  the  most  significant  is  the  limitation  on  the  types  of 
boundary  conditions  that  can  be  described.  The  results  that  have  been  obtained  are  for  a  plate  with 
all  edges  free.  Another  limitation  is  the  complexity  of  the  solution.  Comparing  the  results  obtained 
for  the  response  using  the  3-D  stress  approach  to  the  results  from  the  Mindlin  approach,  when  a  load 
is  applied  on  the  surface  of  a  narrow  thick  plate,  the  results  from  the  two  approaches  are  comparable. 
This  suggests  that  the  approaches  of  Mindlin  and  Mindlin  and  Medick  can  provide  the  required 
results  without  having  to  deal  with  the  full  3-D  stress  equations,  which  are  much  more  complex  to 
solve.  Also,  the  agreement  between  the  two  sets  of  results  verifies  the  approach  that  has  been 
developed  here  to  couple  the  in-plane  and  out-of-plane  waves  with  the  excitation  applied  on  the 
surface  of  the  plate. 

Fluid-Loaded  Plate 

The  influence  of  the  fluid  loading  was  investigated  using  the  same  MPF  technique  as  for  the 
thick  plate  analysis.  The  purpose  of  the  analysis  is  to  investigate  the  influence  of  fluid  loading  on 
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the  power  flow  between  coupled  plates  across  a  discontinuity.  The  analysis  is  performed  for  finite 
plates  which  create  some  severe  complexities  in  the  solution.  Furthermore,  since  the  analysis  of  the 
interaction  through  the  fluid  between  coupled  fluid  loaded  structural  elements  had  not  been  done 
in  the  past,  this  analysis  was  restricted  to  thin  plate  systems.  The  objective  is  to  use  the  results 
obtained  here  as  a  building  block  for  future  work  that  would  include  both  the  in-plane  and  the 
out-of-plane  waves  for  a  thick  plate  structure  or  a  shell  structure.  The  extension  to  thick  plates  with 
both  the  in-plane  and  the  out-of-plane  waves  can  be  accomplished  using  the  developed  approach 
as  a  basis.  The  complexity  in  the  solution  arises  because  of  the  boundaries  of  the  plate,  especially 
if  the  boundaries  are  arbitrary.  To  deal  with  the  arbitrary  boundary  conditions,  the  equations  of 
motion  for  the  plate  and  the  fluid  are  spatial  Fourier  transformed  into  the  wavenumber  domain.  This 
eliminates  the  need  for  assuming  a  known  mode  shape.  Having  performed  the  spatial  transform  a 
1  elation  ship  can  be  established  between  the  scattered  pressure  and  the  motion  of  the  plate.  Since  the 
plate  is  finite,  when  taking  the  spatial  transform  of  the  elastic  part  of  the  equation  of  motion  of  the 
plate,  eight  response  related  parameters  are  obtained  in  the  result  of  the  transform.  From  knowledge 
of  the  boundary  conditions  a  set  of  four  equations  can  be  derived.  The  other  four  equations  are 
obtained  from  the  condition  that  due  to  the  presence  of  the  fluid-loading,  the  response  of  the  plate 
must  always  remain  finite  for  all  frequencies,  even  if  the  structural  damping  is  equal  to  zero.  From 
the  two  sets  of  four  equations,  8  simultaneous  equations  are  set  up  and  solved  to  evaluate  the 
unknown  parameters  in  the  transform  of  the  equation  of  motion. 


In  solving  for  the  scattered  pressure  part  of  the  equation  of  motion,  this  equation  only  holds 
over  the  surface  of  the  plate.  Outside  the  plate  surface,  different  parameters  have  to  be  introduced 
in  the  equation.  For  example  if  the  plate  is  surrounded  by  an  rigid  baffle,  the  stiffness  term  in  the 
equation  of  motion  outside  of  the  boundaries  of  the  plate  should  be  made  infinite.  The  spatial 
transform  of  the  elastic  part  of  the  equation  of  motion  is  thus  a  truncated  transform  creating  the  eight 
parameters  mentioned  in  the  previous  paragraph.  Mathematically  the  same  truncated  transform 
should  be  applied  to  the  scattered  pressure  component  of  the  equation  of  motion,  the  relationship 
between  the  scattered  pressure  and  the  plate  surface  velocity  is  however  based  on  the  full  transform. 
Therefore  the  spatial  transform  of  the  surface  velocity  cannot  be  directly  substituted  for  the 
truncated  transform  of  the  scattered  pressure.  This  difference  between  the  full  and  truncated 
transforms  is  analogous  to  an  edge  effect.  If  the  scattered  pressure  waves  decay  very  quickly  away 
from  the  edge  of  the  plate  on  the  surface  of  the  baffle,  then  the  differences  between  the  full  and 
truncated  transforms  of  the  scattered  pressure  are  insignificant  For  light  fluid-loading  the 
differences  between  the  truncated  transform  and  the  full  transform,  are  not  significant  Physically 
this  implies  that  the  edge  effects  are  insignificant,  which  is  expected.  Furthermore,  for  light  fluid 
loading  the  solution  including  the  truncation  of  the  scattered  pressure  can  be  readily  obtained  by 
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expanding  the  solution  in  a  taylor  series  which  converges.  The  series  has  a  fluid  density  term  in  the 
numerator  which  for  light  fluid  loading  is  very  small  and  thus  the  solution  converges.  For  heavy 
fluid  loading,  as  expected  the  edge  effects  are  more  significant  The  difference  between  the  full  and 
partial  transform  of  the  scattered  cannot  be  neglected.  In  this  case  however,  expanding  the  solution 
in  a  similar  series  as  in  the  light  fluid  loading  case,  the  solution  diverges.  The  fluid  density  term  in 
the  numerator  is  significant  and  the  terms  of  the  expansion  diverge.  A  solution  including  the 
truncation  can  still  be  obtained,  but  it  is  more  complex  mathematically.  The  complete  solution  is 
in  the  form  of  a  Friedholm  integral  of  the  third  kind,  which  can  be  reduced  to  a  Friedholm  integral 
of  the  second  kind.  The  solution  can  then  be  obtained  using  a  projection  method  where  the  first 
approximation  is  the  solution  without  scattering.  This  approach  has  been  successfully  implemented 
and  a  solution  obtained  for  two  end  coupled  plates  and  two  plates  coupled  to  form  an  L -shaped 
configuration.  In  the  case  of  the  L-shaped  plate,  the  diffraction  from  the  outside  comer  is  not 
considered. 

EXPERIMENTAL  ANALYSIS 
Simulation  Results 

In  the  experimental  phase  of  the  work,  simulation  experiments  were  performed  on  a  thick, 
L-shaped,  finite  plate.  The  L-shaped  plate  configuration  was  selected  since  the  analytical  solution 
which  can  be  modified  to  evaluate  the  components  of  structural  intensity  at  each  location  on  the 
plates,  was  already  available.  Using  the  derived  expressions  from  the  analytical  solution,  structural 
intensity  maps  are  generated  for  each  of  the  components  of  the  wave  motion  in  the  thick  finite  plates. 
An  analysis  was  performed  on  the  influence  of  thickness  and  frequency  or  mode  number  on  the 
structural  intensity  maps.  The  results  of  the  analysis  show  that  the  out-of-plane  structural  intensity 
components  decrease  with  increasing  wavenumber  and  thickness  product  (kh),  for  a  fixed  mode 
number,  while  th  t  in-plane  structural  intensity  components  increase  as  the  kh  value  value  increases. 
Less  vibrational  power  is  imparted  to  the  out-of-plane  motion  with  increases  in  the  thickness  of  the 
plates.  This  result  is  consistent  with  some  of  the  results  obtained  in  the  analytical  section.  Also,  as 
would  be  expected,  for  thin  plates  the  components  of  structural  intensity  associated  with  the  in-plane 
motion  are  insignificant  when  compared  to  the  structural  intensity  components  associated  with  the 
out-of  plane  motion.  As  the  thickness  increases,  out-of-plane  motion  structural  intensity 
components  become  the  dominating  components.  The  changes  with  the  mode  number  are  very 
similar  to  the  changes  with  thickness.  The  controlling  parameter  is  the  kh  product,  which  is  modified 
by  either  changes  in  the  thickness  or  the  frequency. 

The  simulation  analysis  also  allowed  for  the  determination  of  the  required  spacing  of  the 
measurement  locations  to  evaluate  the  spatial  transforms  for  the  computation  of  the  structural 
intensity  components  using  finite  differences.  It  is  concluded  that  to  obtain  meaningful  estimates 
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of  the  spatial  derivatives,  spatial  distances  of  the  order  of  1%  of  the  minimum  wavelength  are 
required.  This  can  make  the  separation  distance  very  small.  In  developing  the  simulation  results  a 
separation  distance  for  the  measurements  of  1  mm  on  a  O.S  m  by  1.0  m  plates  is  considered. 

Experimental  Measurements 

Apart  from  the  simulation  experiments,  direct  experimental  measurements  of  structural 
intensity  were  performed  on  a  thick  plexiglass  beam  of  thickness  0.02S  m,  with  one  of  its  ends 
embedded  in  sand,  to  damp  some  of  the  out-of-plane  motion,  and  with  an  external  load  applied  at 
the  free  end.  The  external  load  consists  of  an  in-plane  load  applied  slightly  off  axis  to  induce  both 
in-plane  and  out-of-plane  waves.  The  vibrational  power  associated  with  either  the  in-plane  or  the 
out-of-plane  waves  propagate  at  different  speeds,  with  different  wavenumbers.  Therefore  a 
frequency-wavenumber  decomposition  can  be  used  to  decompose  the  structural  intensity  values 
into  components  associated  with  the  different  wave  types.  The  decomposition  of  the  intensity 
components  is  performed  at  the  cross  spectrum  stage,  before  the  final  intensity  vector  is  computed. 
The  structural  intensity  vectors  can  then  be  plotted  as  a  function  of  frequency  and  wavenumber,  with 
the  magnitude  at  each  frequency  and  wavenumber  value  representing  the  relative  fraction  of  power 
flow  associated  with  the  corresponding  wave  type. 

In  implementing  this  concept  for  the  measurement  of  structural  intensity  on  a  thick  beam, 
vibration  measurements  were  taken  at  a  number  of  equally  spaced  locations  along  the  length  of  the 
beam.  The  in-plane  and  out-of-plane  motions  are  measured  by  a  single  transducer  mounted  on  the 
surface  of  the  beam.  The  in-plane  waves  induce  an  out-of-plane  motion  due  to  the  poisson  effect. 
The  magnitude  of  this  out-of-plane  motion  is  significant  for  thick  structures.  With  the  mainly 
in-plane  excitation  induced  in  the  beam,  the  out-of-plane  motion  associated  with  in-plane  waves 
is  of  the  same  order  of  magnitude  as  the  out-of-plane  motions  associated  with  the  out-of-plane 
waves.  Frequency  cross  spectra  were  formed  for  each  of  the  measurement  locations  and  by 
performing  a  spatial  transform  on  these  measured  cross  spectra,  the  frequency-wavenumber 
structural  intensity  spectrum  is  obtained.  In  the  case  of  the  beam,  in-plane  longitudinal  waves  and 
out-of-plane  transverse  waves  are  the  predominant  wave  motions,  thus  the  frequency-wavenumber 
spectrum  of  the  structural  intensity  shows  two  distinct  set  of  peaks,  one  set  for  the  out-of-plane 
waves,  where  the  relationship  of  wavenumber  with  frequency  at  the  location  of  the  peaks  shows  the 
characteristic  dispersive  curve,  and  the  second  set  of  peak  corresponding  to  in-plane  longitudinal 
waves,  where  the  propagation  speed  is  independent  of  frequency.  If  instead  of  the  beam  a  plate  was 
considered,  then  apart  from  the  longitudinal  in-plane  waves,  in-plane  shear  waves  would  also  have 
been  present 

Additionally,  measurements  were  performed  on  a  semi-infinite  thick  beam  with  and  without 
a  discontinuity.  The  measurements  of  the  beam  response  were  in  this  case  obtained  using  a 
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stereophonic  magnetic  head  mounted  on  a  scanning  table.  Discrete  cross-spectral  density  data 
between  the  response  at  selected  locations  on  the  surface  of  the  beam  and  the  input  load  were  first 
determined  and  then  spatial  Fourier  transformed  to  obtain  a  frequency  -  wavenumber  cross  spectra 
of  the  response.  The  cross  spectra  are  manipulated  to  obtain  the  frequency  -  wavenumber  plots  of 
the  structural  intensity.  A  stereophonic  pickup  was  used  to  simplify  the  acquisition  of  the  data. 

The  experimental  results  obtained  for  the  semi-infinite  beam  show  the  expected  two 
components  of  the  structural  intensity  which  propagate  with  different  wave  speed,  hence  different 
wavenumbers,  in  the  frequency-wavenumber  plot  One  component  is  associated  with  the  in-piane 
waves  ( fast  wave  speed)  and  the  other  component  is  associated  with  the  out-of-plane  waves  (slower 
wave  speed).  These  results  verify  the  technique  that  the  scattering  of  vibration  power  can  be 
determined  experimentally  by  using  a  frequency-wavenumber  decomposition  of  the  structural 
intensity  flow. 

Signal  Processing 

Work  was  also  performed  on  the  signal  processing  approach  required  to  manipulate  the 
acquired  experimental  data  on  a  two  dimensional,  finite  surface  to  deduce  the 
frequency-wavenumber  spectrum  of  the  structural  intensity.  In  the  case  of  a  two  dimensional 
structure,  wavenumber  components  exist  for  each  dimension  of  the  structure.  For  a  structure  with 
both  in-plane  and  out-of-plane  wave  motion,  wavenumber  components  will  also  exist  for  each 
wave  type  in  every  direction.  The  proper  combination  of  the  components  in  one  direction  with  the 
components  in  the  orthogonal  direction  will  result  in  the  wavenumber  values  for  the  selected 
frequency,  and  for  the  particular  wave  type.  Without  a  priori  knowledge,  the  wave  number 
components  cannot  be  properly  combined,  however  this  problem  can  be  worked  around  by  using 
a  two  dimensional  transform.  In  this  case,  the  peaks  in  the  two  dimensional  wavenumber  spectrum 
will  combine  to  represent  the  proper  values  for  the  wavenumbers  for  each  of  the  existing  wave  types 
in  the  structure.  For  example  if  a  wavenumber  spectrum  is  computed  at  a  fixed  frequency,  for  a  finite 
plate  with  both  in-plane  and  out-of-plane  motion,  then  the  two-dimensional  spectrum  will  have 
three  peaks,  one  corresponding  to  the  out-of-plane  wavenumber,  another  for  the  in-plane 
longitudinal  wavenumber  and  the  third  for  the  in-plane  shear  wavenumber.  If  multiple 
two-dimensional  spectra  are  obtained  as  a  function  of  frequency,  the  trace  of  the  peaks  in  these 
spectra  will  represent  the  dispersion  characteristics  of  the  respective  wave  type.  These  result  have 
been  obtained  for  the  structural  intensity  frequency-wavenumber  spectrum  for  a  finite  thick  plate 
structure,  using  simulated  measured  data. 

An  issue  which  is  of  concern  when  performing  spatial  data  processing  with  finite  structures, 
especially  for  low  order  modes,  is  the  proper  identification  of  the  wavenumber  components  given 
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the  sever  truncation  of  the  vibration  motion  shapes.  In  the  work  performed  here,  to  better  distinguish 
wavenumber  components,  the  measured  data  is  extended  by  either  duplicating  the  data  using 
replicas  of  the  u measured  data  or  by  zero  padding.  Padding  with  zeros  increases  the  resolution  but 
does  not  eliminate  the  side  band  effects.  On  the  other  hand,  duplication  by  mirror  imaging  the  data 
such  that  discontinuities  are  not  introduced  eliminates  or  minimizes  the  side  band  effects.  Either  of 
these  two  techniques  are  used  in  the  processing  of  the  simulated  data  to  generate  the 
two-dimensional  frequency  wavenumber  structural  intensity  plots.  Using  this  approach,  it  is  shown 
that  [1 1]  a  good  estimate  of  the  wavenumber  of  the  mode  is  obtained.  This  is  demonstrated  by  the 
consistency  of  the  dispersion  characteristics  obtained  for  a  particular  wave  type  with  what  is 
expected  from  theory.  This  issue  will  not  be  a  problem  for  large  structures  or  for  high  order  modes. 

Non-Contacting  Measurement  Technique 

A  laser  based  approach  has  been  developed  to  perform  measurements  of  structural  intensity 
associated  with  different  wave  components.  In  developing  the  approach,  efforts  were  made  to 
determine  not  only  the  motion  of  the  structure  but  also  the  first,  second  and  third  order  derivatives. 
The  reason  behind  this  requirements  are  twofold.  First,  to  eliminate  errors  that  are  typically 
introduced  in  the  evaluation  of  the  derivative  through  multiplying  the  wavenumber  spectrum  of  the 
response  by  the  wavenumber.  Since  the  wavenumber  spectrum  is  derived  from  measurements, 
estimates  of  high  order  derivatives  based  on  these  measurement  have  a  high  level  of  error  [12]. 
Second,  by  directly  measuring  the  derivatives,  real  time  estimates  of  the  structural  intensity  are 
possible  without  further  post  processing,  if  a  wavenumber  decomposition  is  not  required.  If  more 
than  one  component  of  the  structural  intensity  is  present  and  it  is  required  to  discriminate  the 
components  using  a  wavenumber  decomposition,  then  post  processing  is  required  to  obtain  the 
wavenumber  spectra. 

The  approach  used  to  estimate  the  derivatives  is  based  on  finite  differences.  With  this 
approach  it  is  required  that  the  response  at  multiple  locations  is  measured  simultaneously  by  the  laser 
system.  The  simultaneous  measurements  require  a  high  degree  of  accuracy  especially  with  regards 
to  the  phase  of  the  signal.  The  importance  of  phase  accuracy  in  the  measurement  of  the  structural 
intensity  is  well  documented  [13].  The  determination  of  the  derivatives  using  finite  differences 
would  require  the  same  levri  of  accuracy  as  if  the  structural  intensity  is  being  directly  measured. 
Another  aspect  of  the  finite  difference  approach  for  the  determination  of  the  response  derivatives 
is  that  the  separation  distance  between  the  measurement  locations  must  be  less  than  10%  of  the 
minimum  wavelength  for  which  measurements  are  required.  To  be  able  to  extend  the  measurement 
up  to  as  high  a  frequency  as  possible,  the  finite  difference  spacing  is  reduced.  With  a  laser  based 
approach,  since  the  thickness  of  the  laser  beam  is  very  small,  the  spacing  can  also  be  made  very 
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small.  However  this  creates  a  larger  constraint  on  the  maximum  error  that  is  allowed  in  the 
measurement  of  the  amplitude  and  the  phase  between  the  measuring  channels. 


Fiber  optic  techniques  are  commonly  used  in  laser  vibrometers.  However  if  the  fiber  is 
strained,  such  as  would  be  the  case  if  the  fiber  is  bent,  an  unknown  static  phase  shift  is  introduced 
in  the  signal  passing  through  the  fibre  optic  cable.  If  only  one  measurement  is  being  performed  by 
one  single  fibre,  since  tne  phase  shift  is  static,  this  will  not  influence  the  measurement  However 
when  multiple  point  measurements,  with  different  optical  fibers,  are  performed,  a  non-recoverable 
static  phase  error  is  introduced  in  each  fiber  optic  cable.  This  phase  shift  destroys  the  accuracy 
between  the  channels  and  thus  the  relative  measurements  between  the  channels.  Finally  the  laser 
based  system  had  to  be  capable  of  handling  more  than  one  frequency  component.  That  is  for  the 
frequency  wavenumber  scheme  to  discriminate  between  different  wave  components  of  the  structural 
intensity  an  broad  band  frequency  analysis  is  required. 

With  the  above  requirements  established,  a  compact  laser  based  measuring  set-up  with  all 
the  laser  components  right  up  at  the  measuring  front  end,  to  eliminate  as  much  as  possible  static 
phase  and  magnitude  errors  has  been  developed.  The  set-up  would  be  capable  of  measuring  up  to 
the  third  derivative.  However  significant  difficulties  were  encountered  with  the  set-up  and  to-date 
the  set-up  is  still  non  functional.  Outside  companies  were  asked  to  help  in  this  set-up,  (United 
technologies  Optical  Systems  and  Advanced  Acoustics  Concepts),  supplying  expertise  in  optics  and 
high  frequency  electronic  circuitry.  Instabilities  in  the  electronic  circuitry,  caused  by  high  frequency 
interference  on  the  measurement  signals,  that  conditions  the  optical  signal  before  it  is  passed  to  the 
data  acquisition  system  have  not  been  eliminated. 

The  set-up  of  this  laser  consists  of  a  diode  laser  source  mounted  at  the  head  of  the  multi 
channel  system.  The  single  source  is  then  divided  into  1 3  channels  through  a  series  of  aperture  plates 
and  lenslets.  The  front  end  set  of  the  lenslets  will  focus  the  multiple  beams  on  the  vibrating  target 
which  is  5  to  10  cm  away.  The  total  size  of  the  head  which  includes  all  the  optics  is  0.33  m  by  0.33 
m.  The  selected  set-up  optimizes  phase  accuracy  between  the  channels.  The  back  scattered  light 
from  the  target  is  collected  by  two  sets  of  detectors.  The  output  of  the  detectors  feeds  to  a  processing 
unit  with  a  one  channel  handling  each  pair  of  detectors.  By  combining  the  signals  from  the  pair  of 
detectors,  the  absolute  phase  of  the  channel  can  be  estimated  to  an  accuracy  of  a  few  milli  radians. 
Whether  the  processing  can  be  done  in  real  time  is  dependent  on  how  fast  the  collected  data  can  be 
transferred  from  the  memory  on  the  processing  boards  for  each  pair  of  detectors  to  the  processing 
system,  and  how  fast  the  processing  system  can  manipulate  this  data  to  compute  the  structural 
intensity  results  and  store  these  results  on  some  medium. 
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This  part  of  the  project,  to  set  up  a  laser  based  measurement  approach  to  measure  the  surface 
response  and  its  derivatives  from  which  both  the  in-plane  and  the  out-of  plane  structural  intensity 
components  can  be  simultaneously  measured  in  frequency  and  wavenumber  has  not  been  very 
successful.  The  main  reason  for  this  is  an  underestimate  of  the  time  and  expense  required  to  develop 
an  untested  approach  and  set-up  an  experimental  measurement  technique  by  which  all  components 
of  the  structural  intensity  can  be  measured.  Development  of  this  approach  and  set-up  will  continue 
even  after  the  end  date  of  this  project.  Given  the  effort  that  has  already  been  put  into  the 
development,  and  the  fact  that  now  only  engineering  type  problems  remain,  work  will  continue  to 
be  able  to  perform  some  measurements.  It  has  been  shown  by  using,  other  forms  of  measurements 
of  structural  intensity,  that  the  data  that  will  be  generated  by  this  laser  based  set-up  can  be  used  to 
discriminate  between  the  intensity  components  associated  with  different  wave  types. 
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FUTURE  WORK 

Analytical  and  experimental  tools  have  been  developed  which  can  be  used  to  better 
understand  the  flow  of  vibrational  power  across  structural  discontinuities  in  thick  connected  plate 
like  structures.  The  work  done  thus  far,  while  establishing  the  basis,  is  by  no  means  complete.  Work 
should  be  continued  in  this  area  to  generate  results  for  different  thick  plate  and  discontinuities 
combinations.  A  more  extensive  experimental  program  is  also  required.  Finally  the  techniques 
should  be  extended  to  cylindrical  shell  structures  including  fluid-loading  effects. 
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Power  Flow  Through  The  Junction 
Between  Two  End  Coupled  Plates 


In  previous  work  [1]  equations  were  developed  for  the  power  flow  across  the  junction 
between  two  plates  joined  in  an  L-shaped  configuration.  In  that  case  the  two  plates  were  joined  at 
right  angle  which  coupled  the  in-plane  waves  with  the  out-of-plane  waves  at  the  junction.  The 
out-of-plane  waves  of  the  first  plate  would  couple  with  the  in-plane  waves  of  the  second  plate  at 
the  junction  and  vice  versa.  The  solution  was  developed  in  terms  of  mobility  functions  which  were 
derived  from  Mindlin’s  thick  plate  equations  for  flexural  out-of-plane  motion  [2]  and  thin  plate 
in-plane  equations  of  motion  in  which  only  the  thickness  stretch  motion  was  considered.  An 
improvement  to  this  analysis  uses  Kane  and  Mindlin’s  [3]  in-plane  wave  equations  of  motion. 
Essentially  these  equations  apply  for  the  mid-plane  of  the  plate  with  assumptions  made  regarding 
the  variation  of  the  displacement  as  a  function  of  the  thickness  direction.  All  forms  of  loading,  both 
for  the  in-plane  and  the  out-of-plane  waves  are  describe  with  respect  to  the  mid-plane  of  the  plate. 
The  consequence  of  this  is  that  in-plane  and  out-of-plane  waves  remain  uncoupled.  This  issue  is 
important  when  trying  to  obtain  a  solution  for  end  coupled  plates. 

In  the  case  of  the  end  coupled  plates  (figure  1),  because  the  junction  between  the  two  plates 
is  in  the  same  plane  as  that  of  the  plates,  there  is  no  angle  which  couples  the  in-plane  with  the 
out-of-plane  waves.  That  is,  if  the  plates  are  subjected  to  a  transverse  load  this  would  only  induce 
out-of-plane  motion.  While  if  the  plate  is  subjected  to  an  in-plane  load  it  would  only  induce 
in-plane  motion.  This  is  contrary  to  what  would  be  expected  physically.  Intuitively  it  would  be 
expected  that  if  a  force  is  applied  perpendicular  to  the  surface  of  a  thick  plate  slab,  the  force  will 
induce  both  in-plane  and  out-of-plane  motion. 
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Various  approaches  were  considered  to  work  around  this  problem  [4].  The  literature  was 
surveyed  and  while  solutions  based  on  3-D  stress  analysis  were  found,  these  were  very  restrictive 
in  their  use,  especially  the  type  of  boundary  conditions  that  can  be  considered.  Using  a  variational 
approach  it  is  shown  that  the  Mindlin  [2]  and  Kane  and  Mindlin  [3],  approaches  for  out-of-plane 
and  in-plane  waves  can  be  simultaneously  derived  from  the  3-D  stress  equations,  using  the 
assumption  of  linear  variation  for  the  deformations  and  the  strains  in  the  thickness  direction  of  the 
plate.  Even  though  the  equations  where  derived  from  the  same  3-D  stress  equations,  the  results  are 
two  sets  of  independent  uncoupled  equations,  one  set  describing  the  in-plane  motion  and  another 
set  describing  the  out-of-plane  motion,  with  the  solution  being  in  terms  of  the  deformation  of  the 
mid-plane  of  the  thick  nlate.  The  deformation  at  planes  away  from  the  mid-plane  are  derived  from 
the  assumed  linear  variation  in  the  thickness  direction. 

Therefore,  the  decoupling  between  the  two  wave  forms  is  a  consequence  of  the  linear 
dependency  in  the  thickness  direction.  Furthermore,  the  two  waveforms  remain  uncoupled  even 
when  external  forces  are  introduced.  This  is  because  the  solutio' ,  temal  forces  and  boundary 
conditions  are  all  being  defined  with  respect  to  the  plate  mid-plane.  If  the  forces  and  boundary 
conditions  are  to  be  defined  relative  to  one  of  the  plate  surfaces  and  the  same  type  of  variation 
assumed  over  the  thickness  as  that  for  the  deformation  and  strain,  then  coupling  would  be  obtained 
between  the  in-plane  and  out-of-plane  motions  [4].  This  approach  for  coupling  the  two  wave 
motions  is  here  implemented  together  with  a  refinement  of  the  variation  along  the  thickness 
direction  by  using  a  Legendre  polynomial  expansion  as  used  by  Mindlin  and  Medick  [5]. 


COUPLED  IN-PLANE  AND  OUT-OF-PLANE  WAVES  USING  MINDLIN’S 
APPROACH. 


In  Mindlin  and  Medick’s  approach  [5],  variations  along  the  thickness  direction  are  expanded 
using  a  Legendre  Polynomial.  That  is, 


U(x,y,z) 

-  V(x,y,z)  . 
W(x,y,  z) 

V  . 
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where  =  1,  P^j  =  ^,andP2^j  =  \  ~  ij.ctc.  These  displacements  represented 

by  the  first  few  terms  of  the  expansion  for  all  the  plate  motions,  including  those  for  the  out-of-plane 
waves,  are  graphically  shown  in  figure  (2). 

The  equations  of  motion  ^om  Mindlin  [2]  and  Mindlin  and  Medick  [S]  are  given  by, 
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where  equations  (2)  to  (6)  describe  the  in-plane  motion  of  the  plate  and  equations  (7)  to  (9)  describe 
the  out-of-plane  plate  motion.  The  solution  to  these  equations  of  motion  is  obtained  by  defining 
two  potential  functions  for  the  out-of-plane  displacement  ( 6  and  V'o)  anc*  five  potentials  for  the 
in-plane  displacements  (0lt  02,  03  tpi  and  02  )• 

That  is  for  the  out-of-plane  displacements  ui  and  w\  can  be  described  by  the  two  potentials: . 

dd  ,  #o 

Hi  =  al  +  1F  10- 

1  dy  dx 

which  when  substituted  into  the  equations  of  motion  results  in  the  following  three  equations  of 
equilibrium. 

Tx  v20  +  (*V  -  s "  *)e  -  s  -  ‘wo]  +  =  o  ,2. 

+  (Rd04  -  S-'Y  -  S-'w0]+  =  0  13. 

V2(0  +  w0)  +  S6q4  Wq  =  0  14. 


Eliminating  $  and  y. 


v2(v2  +  6}y0  =  0 
(v2  +  a,2)(v2  +  322)h.„  =  0 
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from  which  alternative  expressions  for  the  out-of-plane  displacements  are  given  by, 
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For  the  in-plane  displacements. 
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which  when  substituted  into  the  equations  of  motion  results  in  the  following  five  equations  of 
equilibrium. 

Vfyi  +  =  0  27. 
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where  for  a  plate  of  thickness  2h 
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In  equations  (27)  to  (31),  v  is  poisson’s  ratio  and  E  and  G  are  respectively  the  modulus  of  elasticity 
and  the  shear  modulus. 

For  a  plate  simply  supported  along  two  edges  (y  =  0  and  a)  and  free  on  the  other  two  edges 
(x  =  0  or  b),  the  solutions  that  satisfy  the  boundary  conditions  are  given  by 
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e2  =  (c2e^  +  D2e~^  sin(~p 

)  33. 

V'o  =  +  D3e~**jcos(^- 

)  34. 

<pt  =  {a^x  +  Bte~^  cos^^^j 

35. 

=  (a4c***  + 

1  36‘ 

ip2  =  (AjC***  + 

|  37. 

A  . 

where  ki ,  (i  =  1,  2,  3)  are  the  out-of-plane  wavenumber  components  and  kt  (i  =  1  to  5)  are  the 
in-plane  wavenumber  components,  respectively  defined  by 
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41. 

For  a  force  applied  on  one  of  the  surfaces  of  the  plate, 

F(x,  y,  z )  =  F(x,  y)  <5(z  -  /i) 

42. 

This  force  can  be  decomposed  using  the  same  Legendre  polynomial  expansion  as  for  the 
deformation  and  the  strains. 

Fix,  y,  z)  =  ^  Pn(^j  Fn(x,  y)  43. 

i*  0 

The  results  of  this  expansion  with  different  number  of  terms  retained  in  the  expansion  is 
shown  in  figure  (3).  As  the  number  of  terms  increases,  the  result  approaches  a  delta  function  at  the 
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surface  of  the  plate  (z/h  =  1) .  This  expansion  gives  the  coupling  between  the  in-plane  and  the 
out-of-plane  waves  is  obtained.  For  a  force  applied  as  shown  in  figure  (4),  retaining  only  the  first 
two  terms  of  me  Legendre  polynomial  expansion, 

Fix ,  y,  z)  =  Fix ,  y )  +  Fix,  y)|  + .  44. 

That  is,  the  force  is  described  by  a  force  of  equal  magnitude  applied  at  the  mid-plane  of  the 
plate  and  two  forces  with  apposite  direction  applied  on  either  surface  of  the  plate  creating  a 
squeezing  load  [5].  The  transverse  force  applied  at  the  plate  mid-plane  only  induces  out-of-plane 
motion,  while  because  of  the  symmetry  about  the  plate  mid-plane  the  squeezing  force  only  induces 
in-plane  motion.  This  same  type  of  approach  was  used  by  Lyamyshev  [6]. 


Figure  3.  Sum  of  contribution  from  elements  of  the  Legendre  Polynomial  expansion  for  the 
surface  applied  load. - :  5  terms; - :  10  terms; . :  20  terms. 

The  first  few  terms  of  the  expansion  match  the  linearly  assumed  variation  in  the  Mindlin  [2] 
and  Kane  and  Mindlin  [3]  solutions.  The  terms  of  the  polynomial  expansion  describe  both  the 
out-of-plane  wave  motion  (W0,  Ut,  Vj)  and  the  in-plane  wave  motion  (W0,  Uq,  U2,  Vq,  Vj), 
where  W,  is  the  displacement  perpendicular  to  the  plate  surface  and  Uj  and  Vt  are  the  displacements 
in  the  plane  of  the  plate. 

To  check  the  validity  of  this  approach,  especially  the  expansion  of  the  surface  force  using 
the  Legendre  polynomials,  the  solution  for  the  response  of  a  narrow  and  long  plate  (equivalent  to 
a  thick  beam)  with  an  edge  transverse  load  applied  on  the  top  surface  of  the  plate  has  been  obtained 
using  the  above  approach  and  compared  to  a  solution  obtained  a  3-D  stress  analysis  similar  to  the 
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one  used  by  Hutchinson  and  Zillmer  [7].  The  Mindlin  and  Mindlin  and  Medick  plate  equations 
(equations  ( 19)  to  (26))  are  used  in  this  solution  for  the  thick  beam  by  letting  n  in  the  sine  and  cosine 
functions  of  equations  (32)  to  (37)  to  be  equal  to  zero. 


x -  - 

in-plane  out-of-plane 


Figure  4.  First  two  terms  in  the  surface  load  expansion. 

3-D  ELASTICITY  SOLUTION 

Considering  a  3-D  slab  as  shown  in  figure  (5)  where  2a  and  2b  are  much  smaller  than  2c 
(thick  long  beam  equivalent  to  a  plate  with  very  narrow  width)  with  free  boundary  conditions  at  its 
ends  and  with  a  transverse  force  applied  at  one  of  the  ends  on  the  surface  of  the  slab. 


Y 


Z 

Figure  5.  Geometry  for  3-D  elasticity  analysis. 


The  3-D  elasticity  solution  starts  from  the  basic  3-D  stress  equations. 
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a  +  G)  v  (V.u)  +  G(v2  u)  =  e(^|) 


45. 


The  solution  to  this  equation  is  obtained  in  terms  of  Helmholtz  displacement  potentials, 

U  =  V(j)  +  V  X  Tj>  46. 


where  u  is  the  displacement  vector  given  by, 

u  =  (ui  +  vj  +  wk) 


47. 


and  9  and  tp  are  scalar  and  vector  potential  functions  respectively.  These  potentials  must  satisfy 
wave  equations  of  the  form. 


48. 


49. 


where  ct  is  the  dilatational  wave  velocity  and  c2  is  the  shear  wave  velocity. 

The  solutions  for  the  two  potentials  in  Cartesian  coordinates  are  given  by. 
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51. 


The  terms  in  the  3  sets  of  brackets  have  8  combinations,  only  four  of  which  are  independent  The 
solution  for  the  displacement  vector  is  thus 
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The  overall  response,  as  given  by  the  above  equation  is  divided  into  four  components  as  indicated 
by  the  roman  numerals. 

Similar  expressions  can  be  derived  for  the  stresses  solutions  which  for  the  direct  stresses 
these  are  given  by, 

Oj  =  Xe  +  2G6J  53. 

where  e  =  du  J dx,  and  for  the  shear  stresses  these  are  given  by. 

^  54. 


where  Yy  =  du  j dXj  +  dllj/ dx y  The  combinations  of  all  the  solutions  are  shown  tabulated  in 
Table  I  which  is  taken  from  [7]  In  this  table  CL  =  (o2X/c\ . 


Table  I: 

Solution  Components  in  Cartesian  Coordinates 
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The  contributions  from  each  of  the  solutions  of  the  vector  and  scalar  potential  functions,  for 
the  four  combinations,  are  arranged  is  such  a  way  that  a  physical  explanation  can  be  attributed  to 
these  contributions. 

For  example,  if  the  upper  terms  in  each  set  of  brackets  are  considered,  u  is  an  odd  function 
of  x,  an  odd  function  of  y,  and  an  even  function  of  z.  v  is  an  even  function  of  x,  an  even  function 
of  y  and  even  function  of  z.,  and  w  is  even  function  of  x,  an  odd  function  of  y  and  an  odd  function 
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of  z.  This  contribution  for  a  structure  with  one  dimension  much  longer  than  the  others,  corresponds 
to  the  even  out-of-plane  beam  modes.  The  other  combinations  would  correspond  to  odd 
out-of-plane  modes,  even  in-plane  modes  and  odd  in-plane  modes. 

For  each  of  the  selected  combination,  the  solutions  for  the  displacements  or  the  stresses  are 
written  in  terms  of  a  linear  combination  of  3  double  series.  For  example,  for  displacement  u, 

NX  NY 

u  =  2  2  Amn  fmn(z)  sin(am  x)  sin(P„  y) 

m  =  1  n  =  1 

NY  NZ 

+  2  2  Bnp  gnp(x)  sin(p„  y)cos(SpZ) 

n=  1  p  =  1 
NZ  NX 

+  22  Cpm  hpm(y)  cos(8p  z)sin(am  x)  55. 

p=  1  m  =  1 

where 

fmn  =  (- am  cos(8z)  +  [-  A2  pn_A3  5]  cos(5z))  56. 

gnp  =  y-  Bt  a  sin(ax)  +  [  -  B2  -  B3  6p  ]  sin(ax>)  57. 

hpm  =  (~  Ci  am  sin(Py)  +  [-  C2  P  -  C3  5p  ]  sin(Py)j  58. 

The  singly  subscripted  constants  and  the  values  of  am,  Sm,  pn,  |5n  8P  and  5P  are  chosen  to 
satisfy  the  boundary  conditions. 

For  free-free  boundaries,  that  is  all  surfaces  free  of  any  stresses,  the  boundary  conditions  arc 
given  by  Tyx(a,y,z)  =  0  ,  Xxz(a,y,  z)  =  0  and  Ox(a,y,  z)  =  0  for  the  surfaces  in  the 
x-plane,  that  is  the  surfaces  perpendicular  to  the  x-axis  (figure  5).  Similar  conditions  can  be  defined 
for  the  surfaces  in  the  y  and  z  planes.  The  6  shear  stress  boundary  conditions  are  satisfied 

term-by-term  by  selecting  the  parameters  am,  am,  p„,  pn  8P  and  $p  to  satisfy  the  conditions. 
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from  which 


a 


62. 

63. 


and  similar  expressions  can  be  derived  for  /?,  ft,  <5,  <F. 


For  the  singly  subscripted  constants,  arbitrary  functions  are  selected  that  must  also  satisfy 
the  6  shear  stress  boundary  conditions.  The  functions  for  the  singly  subscripted  constants  selected 
by  Hutchinson  and  Zillmer  [7]  for  the  case  of  a  plate  or  beam  slab  with  all  edges  free  are  as  follows. 
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The  remaining  3  normal  stress  boundary  conditions  are  satisfied  by  orthogonalization  over 
the  area  of  the  respective  faces  on  which  they  act  If  no  normal  external  load  is  applied  to  a  particular 
surface,  then  the  integral  from  the  orthogonalization  over  the  area  of  the  surface  is  equal  to  zero. 
If  an  external  normal  load  is  applied  on  the  the  surface,  the  orthogonalized  normal  stress  on  the 
loaded  surface  is  equal  to  the  orthogonalization  of  the  external  load. 


For  an  external  normal  load  applied  to  the  surface  in  the  y-plane  ( force  along  the  y-axis 
direction),  the  integral  of  the  stress  function  multiplied  by  the  corresponding  cosine  or  sine 
functions,  over  the  surface  in  the  y-plane  is  equal  to  the  integral  of  the  applied  external  load,  that 
is 

a  c  a  c 


J  J  oy(x,  b,  z)  cos^y)  cos(6jz)  dz  dx  =  j  J  q(x,z)cos(aiX)cos(6jz)  dz  dx  65. 
oo  oo 


where  q(x,z)  is  the  external  load  applied  on  the  surface  y  =  b.  If  the  external  applied  load  is  a  point 
load  of  magnitude  F0  applied  at  (xq,  z0),  then 
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q(x,  z)  =  F06(x  -  Xq)S(z  -  Zq) 


66. 


This  procedure  generates  a  system  of  linear  equations.  The  number  of  equations  generated 
by  the  loading  in  each  of  the  three  orthogonal  directions  is  equal  to  the  number  of  modes  considered 
in  that  directions.  Combining  all  the  linear  equations  a  system  of  (NX*NY  +  NY*NZ  +  NZ*NX) 
equations  is  obtained,  where  NX,  NY  and  NZ  are  the  number  of  modes  in  the  x,  y  and  z  directions 
respectively.  The  solution  to  this  system  of  equations  yields  the  values  for  the  Amn,  Bnp  and  Cpm 
terms  in  equation  (55).  Having  solved  for  these  coefficients,  the  displacements  or  stresses  can  be 
computed  using  equations  of  the  form  of  equation  (55). 

Results  for  the  response  of  a  long  thick  beam  have  been  obtained  using  the  above  approach 
and  the  modified  Mindlin  and  Mindlin  and  Medick  approach.  The  results  are  shown  in  figure  (6). 
There  are  some  discrepancies  in  the  results,  but  these  mainly  stem  from  scaling  factors  and  an 
insufficient  number  of  modes  in  the  three  orthogonal  directions  (see  equation  55).  What  is  of 
importance  is  the  fact  that  in  relative  levels  the  solutions  from  the  two  approaches  match  in  as  far 
as  the  shape  and  ratio  of  in-plane  response  to  out-of-plane  response. 


Figure  6.  Ratio  of  out-of-plane  response  to  in-plane  response. _ :  result  from  approach 

based  on  Mindlin’s  solutions; . :  result  from  approach  based  on  3-D  elasticity  analysis. 

Mindlin  [2j,  and  Mindlin  and  Medick’s  [5]  approach  is  much  more  efficient  than  the  3-D 
stress  analysis  approach  [7]  and  can  be  adapted  to  a  number  of  different  boundary  conditions. 
Therefore  this  method  is  used  to  analyze  the  flow  of  vibrational  power  across  the  junction  between 
two  plates. 
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POWER  FLOW 


To  evaluate  the  power  flow  across  the  junction  between  two  end  coupled  plates,  (figure  1) 
the  solutions  in  equations  (19)  to  (26)  are  used  to  obtain  the  mobility  functions  that  relate  the 
different  motions  (W0,  Wlf  Uq,  Uj,  U2,  Vq,  Vj,  to  the  all  possible  types  of  forces  (T^,  T£,  T^, 
T},  Tj,  Tj,  T5,  Tj)  where  the  same  convention  as  the  one  used  in  [5]  is  used  here  and  extended  to 
include  the  forces  associated  with  the  out-of-plane  motion  (see  figure  7). 


It  should  be  noted  that  the  discontinuity  at  the  junctions  between  the  two  plates  considered 
in  this  analysis  is  a  change  in  the  thickness  between  the  two  plates.  The  discontinuity  does  not 
include  an  external  mobility  (or  impedance)  element  which  when  attached  to  the  surface  of  the  plate 
would  create  a  surface  force  due  to  the  relative  motion  between  the  plate  and  the  discontinuity.  It 
is  important  to  remember  this  point  when  analyzing  the  results  since  the  absence  of  any  external 
forces  at  the  junction  will  keep  the  in-plane  and  out-of-plane  motions  uncoupled  at  the  junction. 

To  simplify  the  displacement  notation,  a  numerical  subscript  representation  will  be  used, 
where  subscripts  1,2  and  3  will  denote  the  x  ,y  and  z  directions  respectively.  That  is  the  response 
in  the  W0  direction  will  be  represented  by  the  subscript  30. 

Since  the  interaction  between  forces  and  displacements  or  velocities  can  only  take  place  if 
both  the  forces  and  the  motions  belong  to  the  same  type  of  wave  motion,  -  in-plane  or  out-of-plane 
-  and  there  are  3  out-of-plane  velocities  (or  displacements)  with  three  corresponding  force 
directions  and  5  in-plane  velocities  (or  displacements)  with  five  corresponding  force  directions, 
then  9  plus  25  mobility  functions  have  to  be  generated.  Furthermore  the  solution  for  the  in-plane 
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and  out-of-plane  waves  can  be  treated  separately.  Therefore,  to  solve  for  the  power  flow  across  the 
junction  between  two  end  coupled  plates,  two  sets  of  matrix  equations  can  be  formulated.  One  set 
solving  for  the  junction  in-plane  motions  and  forces  and  the  other  set  solving  for  the  out-of-plane 
motion  and  forces.  The  two  set  of  matrix  equations,  respectively  for  the  out-of-plane  waves  and 
the  in-plane  waves,  are  as  following, 


^2/3011  +  ^3/30 11 
^2illll  +  ^3illll 
Af  2,21  li  +  ^3/2111 


^2/3051  +  ^3/3051  M2i3060  +  ^3/3060 

T1 

1  l 

“  ^21f3060 

M2iUSl  +  ^3/1151  ^2ill60  +  ^3ill60 

T1 

i5 

= 

~  ^21/1160 

^2i2lSl  +  ^3/2151  ^2i2160  +  ^3/2160 

“  ^21f2160 

*0 

^313011 

^3/3051 

^3/3060 

T1 

1  1 

u, 

= 

^3/1111 

^3/1151 

^3/1160 

A 

V'j 

^3i2111 

^3/2151 

^3:2160 

A 

68. 


M-iauo  ■+■  A/3,3,,0  A/2,31,2  3.31 12  A/ 2,3150  +  A/3,3150  A/23132  A/ 3, 3,52  A/ 2,3161  +  Af 3,3)5] 

7f 

~  A/2i, 3161 

Agnolo  +  A/3,3010  A/2,,012  +  A/3lloi2  A/ 2,)  050  +  A/3,,050  A/^ojj  +  M 3,1052  A/2,,061  +  A#3l|0e, 

7t 

—  A/2|,|06l 

A/2,1210  +  A/  3,32)0  Afj.,212  +  A/3,1212  A/2,1250  +  A/ 3, 3  250  A/211252  +  A/ 3,3  252  Af  2,3261  +  A/3l,2 61 

71 

= 

—  A/21,1261 

A/2.2OIO  +  A/3,2010  A/2,2012  +  A/3,2012  A/ 2,2050  +  A/3,2050  A/ 22052  +  A/  3,2052  A/2,2061  +  A/3,2061 

71 

~  A/212061 

A/2,2210  +  A/3,2210  A/ 2,2212  +  A/3,2212  A/2,2230  +  A/ 3,2252  A/22252  +  MJi2252  A/2,2261  +  A/3,2261 

71 

~  A/2,2261 

L 

69. 


vr 

A/3i3U0 

A/  3,1010 

A/31I21O 

A/  3,2010 

A/  3.2210 

7?" 

t/o 

A/3,3112 

A/ 3, 1012 

A/3, 3212 

A/  3, 2012 

A/3,2212 

71 

t/2 

= 

A/3,3150 

A/3,1050 

A/3,3250 

A/312030 

A/ 3, 2250 

71 

A/3,3152 

A/31IO52 

A/3,3  252 

A/ 3,2052 

A/3.2252 

71 

^2 

A/313I6I 

A/3,1061 

A/3,3261 

A/3,2061 

A/3,2261 

n 

where  the  notation  used  for  the  mobility  terms  consists  of  the  first  number  indicating  the  location 
on  the  structure  for  the  response  and  the  input  load  if  its  an  input  mobility  function  or  the  response 
location  followed  by  the  force  location  if  its  a  transfer  mobility  function,  following  is  an  i  or  a  t  to 
indicate  whether  the  mobility  function  is  an  input  or  transfer  mobility  and  then  the  numerical 
representation  of  the  motion  followed  by  the  subscript  and  superscript  of  the  force.  For  example 
^2i3060  represents  an  input  mobility  at  location  2  for  the  30  response  due  to  a  force  7^  and 

would  represent  the  transfer  mobility  for  the  31  response  at  location  2  due  to  a  force  T\  at  location 
1 .  The  combination  of  the  result  for  the  force  vectors,  from  equations  (67)  and  (69),  and  the  velocity 
response  vectors  from  equations  (68)  and  (70)  will  give  the  transferred  power  across  the  junction 
between  the  two  plates  by  the  in-plane  and  out-of-plane  waves. 
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RESULTS 


The  input  and  transferred  power  flow  results  for  two  different  thickness  ratios  are  shown  in 
figures  (8)  (9)  and  (10).  An  interesting  observation  from  these  results  is  that  while  changes  in  the 
plates  thickness  rado  changes  the  modal  frequency  of  the  out-of-piane  waves,  the  change  in 
thickness  rado  has  no  influence  on  the  modal  frequencies  of  the  in-plane  waves.  The  change  in  the 
thickness  changes  the  mass  and  bending  sdffness  and  thus  influences  the  out-of-plane  motion  of 
the  plates.  However  the  modal  frequencies  for  the  in-plane  waves  are  mainly  controlled  by  the  width 
and  length  of  the  plates  which  is  not  modified  in  the  above  analysis.  The  power  input  and  power 
transfer  associated  with  the  out-of-plane  waves  shows  variations  with  thickness  while  the  power 
input  and  power  transfer  associated  with  the  in-plane  waves  remains  in  general  the  same. 
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Figure  8.  Po\-er  Flow  for  the  coupled  plate  problem  for  a  thickness  ratio  of  05. _ :  Input 

power  to  out-of~plane  waves; _ :  transferred  power  by  out-of-plane  waves; . :  Input 

power  to  in-plane  waves; _ :  transferred  power  by  in-plane  waves. 


Another  observation  which  is  the  result  of  the  previously  mentioned  fact  that  there  is  no 
coupling  between  the  in-plane  and  out-of-plane  waves  can  also  be  seen  from  the  results.  This 
decoupling  between  the  in-plane  and  out-of-plane  waves  is  indicated  by  the  fact  that  while  the 
power  flow  associated  with  the  out-of-plane  waves  type  changes  with  the  thickness  ratio,  that 
associated  with  the  in-plane  waves  type  remains  the  same.  The  in-plane  waves  are  induced  in  the 
plate  by  the  application  of  a  transverse  load  because  of  that  the  load  is  applied  on  the  surface  of  the 
plate. 
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Frequency  Hz 

Figure  9.  Power  Flow  for  a  thickness  ratio  of  1.0. _ :  Input  power  to  out-of-plane  waves; 

_ -  transferred  power  by  out-of-plane  waves; Input  power  to  in-plane  waves; 

...  .  transferred  power  by  in-plane  waves. 
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At  the  junction,  since  the  solution  is  based  on  the  mid-plane  of  the  plates,  the  edge  effects 
due  to  the  unequal  thickness  are  not  taken  into  consideration.  There  are  therefore  no  external  or 
induced  forces  at  the  junction  between  the  two  plates  which  would  couple  the  out-of-plane  waves 
with  the  in-plane  waves.  Out-of-plane  wave  energy  or  power  is  not  scattered  into  in-plane  waves 
energy  or  power  and  vice  versa  because  of  the  presence  of  the  junction  which  has  been  selected  here. 
This  will  not  however  be  the  case  if  the  discontinuity  at  the  junction  was  a  bulkhead  or  stiffener 
which  due  to  its  mobility  will  induce  a  force  on  the  surface  of  the  two  plates.  In  this  case  the 
magnitude  of  this  force  will  be  dependent  on  the  combined  motion  at  the  plate  surface  which  thus 
couples  the  two  wave  types.  If  the  bulkhead  or  stiffener  had  only  lateral  stiffness  in  the  direction 
perpendicular  to  the  plate  surface  and  the  mobility  of  the  bulkhead  is  MB,  then  the  induced  force 
at  the  junction  will  have  a  magnitude  of  (W0  4-  WX)/MB.  Since  this  force  is  applied  on  one  of  the 
plate  surfaces  the  two  waves  would  be  coupled.  In  this  case  a  solution  similar  to  the  one  developed 
in  [8]  can  be  applied. 

CONCLUSIONS 

A  solution  for  the  power  flow  across  the  junction  between  two  plates  has  been  developed 
based  on  Mindlin  and  Mindlin  and  Medick  approaches  for  the  the  response  of  a  thick  plate.  In  these 
approaches  a  Legendre  Polynomial  expansion  is  used  in  the  thickness  direction  to  simplify  the 
governing  equations  of  motion.  Solutions  are  obtained  for  the  equations  of  motion  in  the  mid-plane 
of  the  plate  with  the  Legendre  Polynomials  describing  the  distribution  of  the  displacements  in  the 
thickness  direction.  The  components  of  this  expansion  are  orthogonal  and  the  consequence  is  that 
the  in-plane  and  out-of-plane  wave  equations  become  uncoupled. 

To  couple  the  in-plane  and  out-of  plane  waves  due  to  the  actions  of  an  externally  applied 
load  on  the  plate  surface,  the  load  is  also  expanded  into  components  using  the  same  Legendre 
Polynomial  expansion.  The  results  from  this  approach  are  verified  by  comparison  to  results 
obtained  using  3-D  elasticity  equations  without  simplifying  assumptions.  The  solution  for  the 
out-of-plane  and  the  in-plane  wave  motions  is  developed  in  this  case  based  on  a  superposition  of 
dilatational  and  shear  motions.  General  agreement  is  obtained  for  the  ratio  of  the  out-of-plane 
response  to  the  in-plane  response  between  the  approach  based  on  Mindlin ’s  solutions  and  the  3-D 
elasticity  approach.  There  are  some  discrepancies  in  the  the  agreement  between  the  2  techniques, 
but  this  is  mainly  attributed  to  insufficient  number  of  modes  considered  in  the  3-D  elasticity 
solution.  The  number  of  terms  required  to  obtain  adequate  representation  of  the  physical  behavior 
requires  the  size  of  the  matrix  equation  to  be  very  large,  requiring  significant  computation  time.  The 
approach  based  on  Mindlin ’s  in-plane  and  out-of-plane  wave  solutions  is  computationally  much 
more  efficient 
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When  using  Mindlin’s  approach,  without  external  loads  -  direct  loading  or  loading  from 
externally  attached  structures  (stiffeners  or  bulkheads)  -  the  in-plane  waves  and  out-of-plane 
waves  remain  uncoupled  at  the  junction.  This  is  a  deficiency  of  the  model  due  to  the  fact  that  the 
solution  is  described  relative  to  the  plate  mid-plane.  In  the  case  of  the  two  end  coupled  plates,  with 
different  thickness,  with  external  transverse  load  applied  on  one  of  the  plate  surfaces,  vibrational 
power  is  input  to  each  wave  type,  however  no  scatter  from  one  wave  type  to  the  other  occurs  at  the 
junction.  Coupling  between  the  in-plane  and  out-of-plane  waves  -  which  may  result  in  the 
scattering  of  vibrational  power  from  one  wave  type  to  the  other  -  would  be  obtained  at  the  junction 
if  the  junction  includes  a  discontinuity  which  would  create  loading  on  one  of  the  plate  surfaces. 
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Acoustic  Scattering  From  Fluid -Loaded  Plates 
With  Arbitrary  Boundary  Conditions. 


Acoustic  radiation  from  infinite  fluid  loaded  plates  have  been  extensively  considered  in  the 
past  [1,2].  Similarly  for  finite  fluid-loaded  plates  where  the  boundary  conditions  allow  for  a  modal 
decomposition  of  the  plate  response  [3].  In  other  work  which  deals  with  fluid-loaded  plates  where 
the  boundary  conditions  are  arbitrary,  generally  the  influence  of  the  scattering  is  neglected.  In  this 
paper,  the  dynamic  response  and  acoustic  scattering  from  a  finite,  fluid-loaded,  rectangular  plate, 
simply  supported  on  two  edges  and  with  arbitrary  boundary  conditions  on  the  remaining  edges  and 
set  in  an  infinite  rigid  baffle,  is  considered.  It  is  assumed  that  compressible  fluids  are  present  on  both 
sides  of  the  plate  and  the  plate  is  excited  by  an  obliquely  incident  acoustic  plane  wave.  Because  of 
the  complexity  in  the  solution,  in  this  paper  the  plate  is  considered  a  thin  plate  with  no  in-plane 
waves  and  negligible  shear  and  rotary  inertia  effects.  In  future  work  a  thick  plate  will  be  considered. 
The  work  presented  here  is  a  step  towards  the  analysis  of  the  power  flow  between  two  coupled  plates 
including  both  the  mechanical  and  acoustic  coupling. 

4.  GOVERNING  EQUATIONS 

Consider  a  rectangular  plate  set  in  the  plane  z=0,  bounded  by  0  £  x  <  a  and  0  <  y  <  b, 
and  set  in  an  infinite  rigid  baffle.  The  plate  is  assumed  to  be  simply  supported  along  two  of  the  edges 
(x  =  0  and  a ),  with  arbitrary  boundary  conditions  on  the  remaining  two  edges.  The  semi-infinite 
spaces  z  >  0  and  z  <  0  are  occupied  respectively  by  compressible  fluids  with  density  and 
sound  speeds  q  and  c2.  The  plate  is  characterized  by  a  mass  per  unit  area  D,  thickness  h  and  bending 
stiffness  B.  An  incident  plane  wave  with  angles  of  incidence  <p  and  6  to  the  z  and  x-axes  respectively, 
excites  the  plate  into  flexural  motion  (figure  1). 
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For  a  thin  plate,  the  bending  stiffness  is  given  by  B  =  Eh3/ 12(1  -  v2),  where 
E  =  £0(1  +  jrf)  is  Young’s  Modulus,  rj  is  the  structural  damping,  and  v  is  Poisson’s  ratio. 

The  incident  wave  obeys  the  Helmholtz  equation  in  the  half  space  z  >  0.  The  incident 
acoustic  wave  is  therefore  governed  by  the  equation, 

P  =  71 

where  kx  =  ksintpcosd,  ky  =  fcsin^sinfl,  k2  =  £cos0  and  k  =  oj/cj  is  the  acoustic 
wave  number. 

k2  =  kl  +  k2  +  k\  72. 

All  field  variables  are  considered  to  be  time  harmonic  with  the  convention  e1<ot ,  this  time 
dependency  el<ot  is  however  suppressed  in  the  following  derivations. 

The  equation  governing  small  amplitude  flexural  vibrations  of  a  thin  elastic  plate  in  the 
region  OSjcSci  and  0  S  y  S  b,  is  given  by, 

*W(x,y)  -  D(02W(x,y )  =  [/>]  73. 

where  W(x,y)  is  the  transverse  displacement  of  the  plate  at  position  (x,y) ,  V  4  is  the  square  of  the 
laplacian  operator  and  [P]  is  the  total  pressure  acting  on  the  plate  due  to  the  incident  and  reflected 
acoustic  waves  and  the  scattered  waves, 

[P]  =  [/*(*.?.  2)]z,0_  -  [pfay^z]  +  Pr(x,y,z)  +  Ps^x,y,z)]z^  74. 

P/x,y,z)  and  Pr(x,y,z )  are  respectively  the  incident  and  the  reflected  wave  pressures  if  the  plate 
is  rigid  and  Ph  (x,  y,  z)  and  /*Sj( x,  y,  z)  are  the  scattered  pressure  components  satisfying  the  Helmholtz 
equation  and  the  radiation  condition  [3].  It  has  been  shown  [3]  that  the  reflected  pressure  from  a 
rigid  surface  is  equal  to  the  incident  pressure,  and  therefore  equation  (4)  can  be  simplified  to, 

[^]  =  [Ps2(x>y,z)\Z3S o-  -  [2Pfay,z)  +  />Jl(^>,z)]2==0+  75. 

In  determining  the  response  of  the  plate,  only  the  value  of  the  pressure  on  the  plate  surfaces 
is  considered,  that  is  all  pressure  components  are  evaluated  at  z=0.  and  therefore  for  simplicity  the 
z  dependency  will  be  dropped  from  the  notation  used  in  the  analysis.  At  the  plate-fluid  interface 
the  normal  particle  velocity  of  the  fluid  must  be  equal  to  that  of  the  plate  [3],  that  is  for  0  <  x  s  a 
and  0  Sysi, 


+  jQlfD2W(x,y)  =  0 
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On  the  rigid  baffle  the  normal  velocity  must  vanish  and  the  gradient  of  the  pressure  must  be  equal 
to  zero. 


If  the  plate  is  simply  supported  on  two  of  its  edges,  the  boundary  conditions  along  the  simply 
supported  edges  are  given  by:  VF(0,y)  =  0,  d2W(0,y)/dx2  —  0,  W(a,  y)  —  0  and 

d2W(a,  y)/dx2  —  0.  The  boundary  conditions  on  the  other  two  edges  are  arbitrary.  However  to 
demonstrate  the  approach  and  to  generate  expressions  which  can  be  used  for  the  analysis  of  two  end 
coupled  plate  with  a  stiff  support  at  the  junction,  it  will  be  assumed  that  the  two  remaining  edges 
are  also  simply  supported  together  with  an  edge  moment  at  y=0,  that  is  W(x,  0)  =  0, 
d2W(x,0)/dy2  =  0,  W(x, b)  —  0  and  d2W(x,b)/dy2  =  0  or  —  M(x)/B  depending  on  whether 
the  edge  is  just  simply  supported  or  loaded  with  an  edge  moment. 

The  solution  is  derived  in  the  wavenumber  domain.  This  is  necessary  to  take  into  proper 
consideration  the  influence  of  the  scattered  pressure.  The  forward  and  reverse  spatial  transforms 
are  defined  by  the  Fourier  transform  pair. 


r  «  r  ® 

F(a)  =  I  f{x)eiaxdx  and  f(x)  =  y-  I  F{a)e~iaxda  78. 

J-  00  ^J-O 0 

5.  METHOD  OF  SOLUTION 

In  the  solution  of  the  equation  of  motion,  equation  (3),  a  modal  decomposition  is  assumed 
in  the  direction  perpendicular  to  the  simply  supported  edges,  that  is  along  the  x-axis  direction 

m-  « 

W(x,y)  =  Y  WJy)  sinffi)  79. 

;n  =  0 


Substituting  into  equation  (3), 


d*Wm(y) 

d*y 


-2k, 


m 


\<PWm(y) 

(fiy 


+  k„Wm(y) 


-  Da>2Wm(y) 


-  [Pix,y)] 
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where  km  =  mjt/a,  is  the  modal  wavenumber  along  the  x-axis.  Multiplying  equation  (10)  by 
sin  ( rmx/a\  and  integrating  from  0  to  a. 


+  (*«4  -  tp“)wm(y)  -  ±[-  2P im(y)  -  P,,Jy)  +  81. 

where  kp  -  D(o2/B  is  the  in-vacuo  plate  flexural  wavenumber  and  Pjjy),  Ps^Jy)  and  PtJy)  are 
respectively  defined  by 


*im(y)  =  1 1  Pi(x,y)sin(yffidx  = 


52 J 


=  §[  P«( x,y)sWZTT)dx 

J  o 


,00  =  |  f  /\(*,;y)sin(^)dx 
J  0 


In  the  above  equations  use  has  been  made  of  the  orthogonality  condition 


c 

0 

if  m  *  ri 

Sm(m^x)sm(nf  )dx  -  ■ 

a 

if  m  =  n 

* 0 

2 

h 

In  equation  (12) 


_  (- 1  - 1 
•'m  — 


—  jr ^ 
'vm  ** 
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83. 
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When  the  acoustic  wavenumber  in  the  x  direction  is  equal  to  the  plate  mode  wavenumber 
in  the  x  direction,  that  is  km  =  ±  kx ,  the  coefficient  Jm  is  not  singular,  but  is  given  by. 


Jm  — 


—  if  k  =  k 

2  11  Kx 

—  —  if  k  —  —  k. 

2  u  "’M 
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Because  of  the  nonlinear  form  of  equation  (11),  there  is  no  exact  simple  analytical  solution 
to  this  equation.  If  the  scattering  pressure  components  can  be  neglected,  then  an  analytical  solution 
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is  possible.  Otherwise  a  numerical  approximation  method  has  to  be  used  for  the  solution.  The 
following  sections  present  solutions  to  the  coupled  equations  of  motion  of  the  finite  fluid-loaded 
plate. 

6.  NEGLIGIBLE  SCATTERING  COMPONENT 

Neglecting  the  scattered  pressure  component,  equation  (11)  can  be  solved  analytically  and 
the  solution  can  be  written  in  the  form. 


WQm(y)  =  Wp(y )  +  Wk(y) 


where  Wp  is  the  particular  solution  of  equation  (11), 


Wp(y) 


~  4/»0  •  Jm 

a  J  b 


eJk>y 

[ki  +  qf  -  k* 


and  Wg  is  the  general  solution  of  the  homogeneous  equation. 


88. 
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Wh(y)  =  Y  xe>K'y  +  Y2e  ~*K'y  +  Y^K*  +  YAe~^  90. 

In  equation  (20),  T,  (  i  -  1  to  4  )  are  arbitrary  constants  which  are  determined  from  the 
boundary  conditions.  K  \  and  K  2  are  the  roots  of  the  dispersion  equation  and  are  given  by 

*1  =  ±  jj*}  +  &  and  K2  =  ±  Jk2p  -  k2n.  KiandK2  are  selected  such  that  the  imaginary  part 
is  positive. 

In  the  case  of  a  simply  supported  plate,  the  coefficients  T,  ( 1  =  1  to  4  )  arc  given  by. 


k]  -  K\ 
4*2/sin(tf1&) 


^24  =  ~ 


k]-K\ 

jsin(Kib) 


(eiW  -  eiK'b)  Wp 


r“  =  "  S  ^  ~  e'JK^  Wp<0) 


Y“  -  4^)  -  **) 
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where  the  additional  subscript  A  is  introduced  for  the  simply  supported  case. 

For  a  plate  with  an  edge  moment  excitation  M(x)  at  the  edge  y=0,  the  coefficients 
Yt  ( i  —  1  to  4  )  are  given  by, 

Kg-j^b  KejK,b 

y  _  B  y  _  B 

m  4k2  jsin(Kxb) '  m  4k2jsin(Klb)' 

!±e-)Ki>  &ejKJ> 

4*2  j  s'm(K2b)  ^  4 k2  j  s'm(K2b) 

the  subscript  M  is  to  indicate  the  presence  of  the  edge  moment  In  equation  (22), 


Other  boundary  conditions  can  be  considered  using  a  similar  approach. 

7.  SOLUTION  INCLUDING  SCATTERING 

When  the  scattered  pressure  components  are  not  neglected,  a  relation  between  the  scattered 
pressures  and  the  plate  displacement  is  required  in  the  solution.  From  equation  (6),  the  relation 
between  the  displacement  W  and  the  scattered  pressure  in  the  wavenumber  domain  is  given  by 

Y’PAa'.P)  =  je<D2W(fii',P')  94. 

where  y'  is  the  acoustic  wavenumber  in  the  z-axis  direction 

y'  =  ±  Jk2  -  a'2  -  f}'2  95. 

The  appropriate  sign  and  phase  are  determined  by  the  physics  of  the  problem.  Fourier  transforming 
equation  (11) 

(|82  +  i kjf  -  -  {P2  +  2 km2)[wm\by>*  -  Wm'(0)}  + 

jp(p2+2km2fwm(b^  -  HUO)]  -  jf>[wm"(b)e‘*  -  »V'(0)]  + 
W„'"(b)ei *  -  Wm'"( 0)  =  i[  -  2 P,Jfi)  -  PSJJ3)  +  P,J0)]  96. 
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where 


The  W  primed  terms  in  equation  (26)  are  a  result  of  the  fact  that  the  plate  is  finite.  Introducing  the 
following  simplifying  notation. 


A(fi)  =  (/S2  +  k2mf  -  4 
Ax</J)  -  {p2  +  2ki)e^t 
A20)  =  -  (/32  +  2J&) 


A3<0)  =  -  e^b  A6(fi)  =  ;i8e/W 

A4(J5)  =  1.  An(fi)  =  -  jfiifi2  +  2k2]e^b 

A,(fi)  =  -  jP  At(fi)  =  y/3(/32  +  2k2m) 


98. 


and 


X,  =  wm‘(b) 
X2  =  Wm'(0) 
X3  =  Wm'"(b) 


X 4  =  WV"(0) 

X,  =  w„”(0) 

x6  = 


X7  =  Wm(b) 
Xg  =  wm(0) 


Substituting  into  equation  (26), 

A<j3)Wm<y3)  =  Q(fi)  +  -  A,„<0)] 


100. 


where  the  function  Q(fi)  is  defined  by. 


1  =  8 

m  =  - 1  P/.0)  +  ioi- 

;  =  i 

From  equation  (27)  the  Fourier  transform  of  the  pressure  components  also  ranges  from  0  to 
b.  Therefore,  the  relationship  is  not  in  the  form  shown  in  equation  (24).  The  relationship  between 

P0)  and  Wm(fi)  (derived  in  Appendix  A),  is  given  by, 


Ai 


/■  +  ® 


102. 
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—  p 
B  rs'~ 


03) 


r  +  ® 

=  -a2 

J  -  00 


^'X?03  -  PyWJfi'W 


103. 


where  =  jQ1(o2/(2jt2aB),  X2  =  jQtfi>2/(2n2aB),  G(0)  =  -  1/00)  and  71<t  and  72.  are 

defined  in  Appendix  A. 

Due  to  the  finiteness  of  the  plate,  for  each  wavenumber  0,  the  scattered  pressure  component 
is  dependent  on  all  of  the  plate  displacement  components.  The  fluid-plate  coupling  is  represented 
by  the  parameters  Ai  and  Xj  and  the  functions  Jlm  and  J2m-  The  coefficients  Xi  and  X2  are  functions 
of  the  fluid  loading  parameters  [l]ai  and  a2  which  are  given  by,  a  j  =  o)2q1/B  and  n2  =  o)2q2JB. 

The  functions  Jim  and  represent  the  importance  of  the  contribution  of  each  plate 
wavenumber  to  the  scattered  pressure.  In  particular,  the  closer  the  wave  number  is  to  the 
wavenumber  km,  the  more  significant  is  the  contribution,  (see  Appendix  A). 

Substituting  equation  (32)  and  (33)  into  equation  (30), 

=  Q(fi)  +  K(W„m  104. 


where  K  is  defined  by, 

K(<P)(P)  =  -  [  [Vi Jfi')  +  V2„O3')]g(0  -  pwp'W  105. 
J  —  OO 

The  integral  equation  (34)  represents  the  complete  solution  to  the  plate  scattering  problem. 
In  evaluating  this  equation,  it  is  still  necessary  to  obtain  the  eight  unknown  functions 

Xj  (i  =  l  to  8  ),  which  are  present  in  the  function  Q{p).  Eight  simultaneous  equations  are 
required  to  solve  for  these  unknowns.  Four  of  the  simultaneous  equations  can  be  obtained  from  the 
condition  that  whenever  the  polynomial  A(fi)  goes  to  zero,  the  first  term  in  equation  (34)  also  goes 

to  zero  and  therefore  the  second  term  in  equation  (34)  must  also  vanishes.  Since  4(8)  is  a  polynomial 
of  degree  four  this  condition  gives  four  linear  equations.  The  other  four  equations  come  from  the 
boundary  conditions. 

Equation  (34)  is  a  Fredholm  equation  of  the  third  kind.  The  main  difficulty  of  solving  this 
equation  lies  in  finding  an  accurate  representation  for  the  operator  K.  Different  techniques  are 
available  to  evaluate  K  [4].  These  include,  (a)  Iterative  methods,  (b)  Quadrature  methods,  (c)  Kernel 
approximation  methods  and  (d)  Projection  methods. 
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The  iterative  methods  are  not  applicable  because  of  the  large  values  of  the  parameters  Xj  and 
X2  The  iterative  algorithm  in  general  does  not  converge.  For  the  quadrature  method  a  large  number 
of  points  is  required  to  cover  the  full  range  of  integration  which  makes  this  technique  very 
computationally  intensive.  Thus,  only  the  Kernel  approximation  method  and  the  projection  method 
may  be  possible  to  use  for  obtaining  a  solution. 

8.  KERNEL  APPROXIMATION  METHOD 

When  the  acoustic  wavelength  is  less  than  the  plate  dimensions,  the  behavior  of  the  plate  can 
be  assumed  to  be  similar  to  that  of  an  infinite  plate.  That  is,  the  function  G  approaches  a  dirac  delta 
function,  G(fi)  =  2ji60).  An  approximation  for  the  operator  K  can  therefore  be  obtained  as, 

K(<PX0)  =  ~  2tf[Vu<0>  +  V2m<0)]<W  106. 

Substituting  equation  (36)  into  equation  (34) 

Wm0)  =  - — ^ - 77  107. 

[a<P)  +  2*(V,„03)  +  V2.<P))j 

When  determining  the  unknowns  in  the  function  four  of  the  unknowns  are  determined 

from  the  condition  that  whenever  the  denominator  in  the  equation  (37)  is  equal  to  zero,  the  numerator 
must  also  vanish.  This  denominator  is  a  sum  of  the  polynomial  A(fi)  and  a  perturbation  function. 
It  can  be  assumed  that  the  zeros  of  this  denominator  will  be  in  the  vicinity  of  the  zeros  Kj  and  K2 
of  the  polynomial  A(0).  The  actual  zeros  can  be  found  by  an  iterative  process. 

9.  PROJECTION  METHOD 

The  Fredholm  integral  equation  of  the  third  kind  (equation  34)  can  be  simplified  to  an 
integral  equation  of  the  second  kind.  This  is  done  by  determining  the  four  unknowns  of  the  function 
Q  which  are  not  determined  by  the  boundary  conditions. 

Considering  the  case  of  the  simply  supported  plate,  from  the  boundary  conditions, 
X5  =  X6  =  X7  =  Xg  =  0  and  equation  (31)  becomes, 

1=4 

QM  =  - 1  +  X*A<£)  108. 

i  =  1 

The  coefficients  Xi  (  1  =  1  to  4  )  can  be  obtained  from  equation  (34),  since  when  (3  is 
equal  to  one  of  the  zeros  of  the  polynomial,  A  (  ±  Kx ,  ±  K2),  the  left  hand  side  of  equation  (34) 
goes  to  zero,  the  right  hand  side  must  also  vanishes.  Thus  a  set  of  four  equations  is  obtained  which 
can  be  solved  for  X,  (  i  =  1  to  4  ),  that  is 
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109. 


i=4 

i-1 

where  JTcan  take  values  ±  Xj  or  ±  K2  ■  Solving  this  system  of  equations  analytically,  the  result 
is  of  the  form 

*1  =  XIA  +  *IS  *3  =  X3A  +  *3S 

X2  =  xu  +  *25  x4  =  XAA  +  X4S 

where  Xia,  X2A.  X3A,  X4A  represent  the  solutions  of  equation  (39)  when  the  scattering  effects  are 
neglected.  These  terms  are  given  by, 

Xu  =  ^{[P,_(K,>  -  PJ-  X,)]sin(X2«  -  [#JK2)  -  P,S- X2)]sm(Klb)) 

111. 

Xu  =  ^{[p,  « :i)e-‘K*  -  P,( -  K,)^b]sm(K2b)  - 

[?im(K2)e-iK*  -  PJ-  ^e^sin^fe)}  112. 

=  2b\{P‘-(K i)  '  +  - 

[p,„(X2)  -  PJ  ~  /C2)](AT?  +  2 kl)  sin(X,6)}  113. 

Xu  =  Jg{[PiSKl)'~iK'b  ~  P<S~  ^K'b][Kl  +  2^)sin(/C2 6)  - 

[pjX2)e*«’  -  PJ-  K2)^](A:;  +  2^)sin(/C,f)))  114. 

Xis,  X2s,  X3S,  X4S  are  the  contributions  from  the  scattering  given  by 

XIS  =  -^±{[X(W'„)(K1)  -  K(W„X -  A'1)]sin(JC2*)  - 

[*(tf JflCj)  -  K(Wm)(  -  X2)]  sin(K,«J  115. 

Xjs  =  ^L{[K(Wm)(Kl)e-^b  -  K(Wm)( -  sin(X2«  - 

[x(H'„)(X2)e -  HWM-  X2)^]sin(X,fc)}  116. 
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=  J7i{[K<"'»XA:l>  ~  Ki)]sin(K2«(*2  +  2*5.)  - 

[K(W„XK2)  -  K(WJ( -  K^sMKMK]  +  2*5,)} 


117. 


X«  =  -  K(W„X-  Kl)e’K‘,’}sm(K2b)(Kl  +  2*5,)  - 

[Ki.W^K^e-i^  -  K(Wm)(  -  ATj)^-4]  sin(JC,6)(Ar5  +  2*5,)}  118. 

where  A  —  —  Ajkp  sin(Kxb)  sin (K2b).  Substituting  these  equations  into  equation  (38) 

4  4 

203)  =  - 1  P,Jfi)  +  X^,03)  +  2*0*1 09)  ns 

i-1  i=l 

If  the  scattering  terms  are  neglected  in  equation  (49),  then 

4 

A(fi)wAm<$)  =  - 1 4o3)  +  2x»m  i2°- 

i=i 

where  WAJfi)  is  the  Fourier  transform  of  the  plate  displacement  when  the  scattered  pressure  is 

neglected,  given  by  the  Fourier  transform  of  equation  (18).  Substituting  equations  (41)  to  (48)  into 
equation  (49),  the  following  expression  (see  Appendix  B)  is  obtained. 


20?)  =  A(fi)WAJfi)  + 


(+  « 

[Vi„0’)  +  y2Jfi')]wm0’)K<fi,P')dP’ 


121. 


where 

Kifi,?)  =  jj-(032  -  K2)  sin(K2b) 

[g^  -  P’)(e~jK'b  -  e>Pb)  -  G(-  Kx-  P')[e>K'b  -  ^b)] 

-  (fi2-  A^sm^^G^  -  -  e>Pb)  -  G(-  K2  -  P'){e>KJ>  -  ^)]} 

122. 

And  substituting  into  equation  (34),  the  Fourier  transform  of  the  plate  displacement  is  then  given 
by  a  Fredholm  integral  equation  of  the  second  kind, 


45 


123. 


Wm<0)  =  WJff)  - 

r  +  « 

■jL  [Vi.<0')  +  V2.(0')][G</9  -  /S')  +  K(fi,0')]WJ0')dp' 

J  —  00 

The  displacement  of  the  plate  then  consists  of  a  displacement  contribution  obtained  by 
neglecting  the  scattering  terms  and  a  corrective  component  introduced  to  account  for  the  presence 
and  influence  of  the  fluid  loading. 

If  the  plate  is  under  an  edge  moment  excitation,  a  similar  procedure  as  above  can  be  followed 
to  derive  the  dynamic  response  of  the  plate.  The  only  change  in  the  result  being  in  the  solution 
WAJfi)  to  the  plate  displacement  when  the  scattered  terms  are  neglected.  This  solution  will  be 
referred  to  as  One  can  show  that  the  constants  Xj,  X2,  X3,  X4  in  are  given  by 

Xw  =  -  |f  [tf ,  sin (K2b)  -  K2  sin(^,6)|  124. 

Xm  =  -^[Ar1cos(AT16)sin<A:26)-A:2cos(A:26)sin(A:16)]  125. 

xm  =  -  If  [at, (k\  +  2k2„)S\n(K2b)  -  K2(K\  +  2*J,)sin(Ar,*)]  126. 

Xiu  =  -  |f  [Ki(*1  +  2k2m)sm(K2b)cos(K1b) 

-  K2(K\  +  aysin^^cos^)]  127. 

where  M  is  the  edge  moment,  given  by  equation  (23)  and  the  dynamic  response  of  the  plate  is  given 
by  an  equation  similar  to  equation  (53)  with  W0Jfi)  substituted  for  WAJfi),  where  W0Jfi)  is  the 
plate  response  when  the  scattered  pressure  is  neglected  and  which  is  equal  to  WAJfi)  in  the  case  of 
acoustic  excitation,  in  the  case  of  edge  moment  excitation,  and  WAJfi)  +  in  the 

case  of  both  acoustic  and  edge  moment  excitation.  W0  Jfi)  can  be  derived  from  equations  (18). 

Equation  (53)  is  the  general  solution  to  the  problem  of  scattering  from  a  plate  using  a 
projection  method.  Projection  methods  include  the  Polynomial  Collocation  method,  Galerkin’s 
method,  the  Least  Squares  method,  and  others  [4].  The  polynomial  collocation  consists  of 
approximating  the  solution  to  the  integral  equation  using  a  set  of  polynomial  basis  and  a  set  of 
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collocation  points.  The  performance  of  this  method  depends  on  the  choice  of  this  set  of  points.  For 
the  problem  considered  here,  the  solution  is  highly  oscillatory  and  the  range  of  integration  is  between 
-°o  and  -h»  ,  making  this  approach  inappropriate.  The  Least  Squares  method  is  also  not  suitable  since 
it  requires  the  calculation  of  multiple  integrals.  Galerkin’s  method  consists  of  decomposing  the 
solution  using  a  basis  and  approximating  the  solution  by  truncating  the  decomposition.  The 
performance  of  the  method  depends  only  on  the  choice  of  the  basis.  This  method  is  the  one  used 
in  the  solution. 

10.  GALERKIN’S  METHOD  AND  CHOICE  OF  BASIS 

If  the  scattered  pressure  is  neglected,  the  solution  for  the  plate  displacement  is  of  the  form 
e*y  (equations  19  and  20),  where  K  can  take  any  of  the  complex  numbers,  ±  Kj ,  ±  K2 ,  Ky.  The 
form  of  the  solution  for  the  plate  displacement  including  the  influence  of  the  fluid-loading  will  in 
general  not  significantly  deviate  from  this  form.  Therefore,  it  can  be  assumed  that  the  solution  will 

still  be  of  the  form  e*y ,  but  with  different  values  of  X.  In  the  wavenumber  domain,  the  solution 
will  have  a  component  of  the  form  G(X  +  /?).  Using  a  Taylor  expansion,  G{K  4-  p)  can  be  written 
as. 


G<X  +  fi)  =  ^ 

n-0 


rdnG(0) 

n\  dnfl 


128. 


Let  <pn(fi)  =  dnG(fi)/dn(3,  these  functions  are  a  basis  of  the  space  of  oscillatory  functions. 
That  is,  any  oscillatory  function  can  be  decomposed  uniquely  in  terms  of  these  functions.  The 
evaluation  of  these  functions  is  presented  in  Appendix  C.  The  plate  displacement  in  the  wavenumber 
domain  can  therefore  be  decomposed  in  the  form 


wm(fi)  =  £  129. 

n  =  0 

where  the  coefficients  (  a„  ,  n  =  0  to  »  )  are  to  be  determined  using  equation  (53).  In  the 

Galerkin’r  method  Wm  is  approximated  using  a  finite  summation,  that  is  equation  (58)  is  truncated 
for  some  value  of  N,  which  is  chosen  based  on  the  desired  level  of  accuracy.  Substituting  equation 
(58)  into  equation  (53) 


N 

W0 ,03)  =  £ 


n  =  0 


130. 


where 
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131. 


1 

m 


+  V2.0')][G</3  -  P')  +  K<fi,P')}<p'(fi')dp' 


The  coefficients  an  (n  =  0  to  N)  can  be  obtained  by  solving  a  set  of  equations  of  the  form 
N 

2*«[<  <pn,<f>l  >  +  <  ^*.0/  >]  =  <  W0m><P,  >  132. 

for  0  ^  /  S  N,  where 


f  +  00 

<  M,  v  >  =  I  u(fi)v  *  (fiyip  133. 

J  —  00 

The  asterisk  indicates  a  complex  conjugate.  Substituting  for  <  <pn,<pi  >,  <  ^n.0/  >  and 

<  W0m,<pl  >  ( see  appendix  C) 

<  >  =  2 »«-  134. 

<  V,n»0/  >  and  <  WQm,<t>t  >  are  computed  using  a  numerical  integration.  When  the 
fluid-loading  effects  are  neglected,  the  term  <  ipn^i  >  in  equation  (62)  can  be  neglected.  The 
solution  to  this  system  gives  an  approximation  to  the  exact  solution  W0m  which  is  given  by  equation 
(50)  or  by  equation  (18).  It  has  been  found  that  N=5  gives  an  excellent  approximation  to  VV0^.  It 
is  therefore  assumed  that  N=5  will  give  a  good  approximation  to  the  exact  solution  when  the 
scattering  terms  are  included  in  the  solution. 

11.  INPUT  POWER 

The  power  input  to  the  plate  is  balanced  by  its  dissipation  through  mechanical  losses  and 
acoustic  radiation.  From  the  equation  of  motion  of  the  fluid-loaded  plate, 

-  2P,(x,y)  =  j|[W(*,y)  -  k*V(x,y)]  +  Ps>(x,y)  -  />„(*, y)  135. 

the  power  flow  into  the  plate  from  the  incident  pressure  is  obtained  by  multiplying  both  sides  of  the 
above  equation  by  V  *  (x,  y)/2  and  taking  the  integral  over  the  plate  surface.  The  total  input  power 
to  the  plate  is  thus  given  by, 
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fl,=  -  ±Rejj  |  2  P,(x,y)V*(x,y)dxdy  ■ 


136. 


The  factor  of  2  in  front  of  the  incident  pressure  is  included  to  take  into  account  the  reflected  pressure 
contribution.  Substituting  for  the  displacement  using  equation  (9)  and  using  the  notation  of  equation 
(12) 


137. 


The  spatial  integrals  in  the  expression  for  the  input  power  can  be  replaced  by  integral  in  the 
wavenumber  domain  using  the  complex  Fourier  spatial  transforms.  Therefore 


1=  OD 

I- 


138. 


Pim  is  the  full-range  Fourier  transform  of  Pjjy)  •  From  equation  (12) 

P,Jy>  =  jp  ^  i»- 

therefore 

Jm  <5 (/S'  +  ky)  140. 

and 

m  =  oo 

Tli  =  -  P0  Irnag[jmV*m(-  141. 

m=  1 

12.  SCATTERED  POWER 

The  scattered  power  is  the  power  transferred  to  the  fluid  from  the  plate  and  is  given  by 


Psfay)  -  PSi(x,y)]v  *  (x,y)dxdy> 


142. 
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In  terms  of  modal  contributions  IJS  can  be  written  as 


ns  = 


143. 


From  Appendix  A,  the  relationships  between  the  full-range  Fourier  transform  of  the  scattered 
pressure  components  and  the  plate  surface  velocity  are  given  by, 

144. 


PsJfi')  -  -  145. 

Fourier  transforming  equation  (73), 


m*t  aa  4  +  « 

n,  =  £&  Re[e,i,.03>  +  e?/2jj3)Jiv„(/3)i2dp 

«=  1  j  -  QO 


146. 


From  Appendix  A,  the  real  parts  of  J\Jfi)  and  J2Jfi)  are  not  equal  to  zero  only  for  the  wavenumber 
satisfying  (01  S  k.  Therefore  only  the  components  of  the  plate  displacement  corresponding  to  these 
wavenumbers  contribute  to  the  radiated  sound  power.  The  range  of  integration  in  equation  (76)  can 
be  reduced  to  -  k  <,  0  <>  k. 


ns  =  2  £2  [  Re[e,/,.03)  +  147. 

m- 1  J -k 

To  be  consistent  with  the  fundamental  law  of  the  conservation  of  energy,  this  scattered  power  must 
always  be  less  than  the  input  power  to  the  plate.  According  to  equations  (74)  and  (75),  the  terms 
that  mainly  represents  the  influence  of  the  fluid  loading  on  the  plate  displacement  are  the  integrals 
J lm(fi)  and  It  is  shown  in  Appendix  B  that  at  a  particular  frequency  of  vibration,  the  form 

and  associated  physical  nature,  of  the  fluid  loading  represented  by  J 1  Jfi)  depends  upon  how  the 
plate  wavenumber  fi  in  the  Y  direction  is  placed  with  respect  to  the  acoustic  wavenumber  k.  For 
10!  2:  k,  the  integral  Jijfi)  is  purely  imaginary  and  is  given  by  equation  (A.  12).  This  imaginary  part 
represents  the  effective  mass  loading  on  the  plate  by  the  reactive  component  of  the  fluid  loading. 
For  |0|  5  k,  the  component  J \Jfi)  has  both  real  and  imaginary  parts.  The  real  part  contributes  to 
the  far  Held  radiation  and  it  represents  the  resistive  component  of  the  fluid  loading.  This  real  part 
is  present  at  all  frequencies  for  a  finite  plate. 
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From  Appendix  A,  equation  (A.  15),  for  any  given  wavenumber  /J,  (  y?|  s  k ),  the  maximum  value 
of  the  real  part  of  J i  Jfi)  is  reached  when  the  mode  wavenumber  km  is  less  than  the  value  of  V,  given 

by  V  =  Jk2  -  f32.  In  this  case,  the  range  of  integration  in  the  integral  giving  this  real  pan  includes 
the  maximum  value  of  the  integrand.  Thus  when  the  wavenumbers  of  the  plate  in  both  the  X  and 
Y  directions  are  less  than  the  acoustic  wavenumber,  that  is  above  critical  frequency,  the  acoustic 
radiation  will  reach  a  maximum. 

Of  interest  here  is  the  radiation  from  a  plate  where  none  of  the  resonant  modes  excited  in  the  plate 
are  radiating  as  efficient  acoustic  radiators.  That  is  when  the  plate  is  excited  below  the  critical 
frequency.  For  a  6  mm  thick  steel  plate  in  air  the  critical  frequency  is  around  2  kHz.  Therefore  in 
the  results  the  frequency  range  is  limited  to  1  kHz. 

13.  MECHANICAL  LOSSES 

The  mechanical  losses  component  of  the  power  flow  includes  both  the  power  lost  due  to  the  internal 
structural  damping  and  the  power  transmitted  across  the  junction  if  the  plate  is  a  component  of  a 
global  structure.  An  explicit  expression  for  the  mechanical  losses  is  given  by. 


n 


M 


00  •  n  cb  r 

=  Y  -jaB 

Z.  Aw 
m= 1  Jo  L 


dAV 


m' 


dy 


®*y  -  2km^dy 


dy1 


[C  -  kp4)v(y)  Vm 


( y)dy 
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Integrating  by  parts 


149. 


where  X, ,  i  =  1  to  8  are  defined  by  equations  (29).  In  the  above  equation,  the  spatial  integrals  can 
be  transformed  into  the  wavenumber  domain  as  follows. 


Wiytfdy 


■m: 


\V(fi)\  2d0 


and  the  Fourier  Transform  of  the  derivative  of  the  velocity  is. 


V'(fi)  =  Jw(x1e*s*  -  X2)  -  jpvtfi) 


Thus 


dViy) 

dy 


p§r)  dy  =  2kj  I  Hx^  ~  x*)  -  ^n|*4» 


Similarly  for  the  second  derivative  of  the  velocity. 


(PVWfePVjy) 
dy2  \  (f-y 


j  dy  = 


r +  ao  2 

^  h(x3^  -  x4)  +  ofiXtt*  -  x2)  -  p2m\  dfi 

J  —  00 


The  mechanically  dissipated  power  is  thus, 

nM  =  £  ^^[XjX2  *  -  x^x,  *]  + 


-  <4-  00 

j  [l/»(x,«*»  -  Xj)  -  /32V0)|2  +  (kj  +  2 kjp  -  */) |v(y})|2  48 


In  equation  (84),  the  terms  with  X5  and  X6  represent  the  transmitted  power  if  the  plate  is  connected 
to  another  plate.  The  rest  of  the  terms  represents  the  dissipation  due  to  the  structural  damping. 

In  the  case  of  a  plate  simply  supported  on  all  four  edges,  equation  (84)  reduces  to. 


m=  1 


j  j^fx,^  -  X2)  -  fi2V(fi)\2  + 
(*m4  -  */  +  2*„2/32)|V0S)|2j49 
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For  a  plate  with  an  edge  moment  at  y  =0  and  with  simple  supports  on  the  remaining  edges, 
where  the  edge  moment  is  due  to  the  influence  of  another  structure  attached  to  the  y=0  edge  of  the 
plate,  the  mechanical  dissipation  is 


nu  =  "X  *a*  + 

mm  1  J 

3^  j  -  xl)  -  +  (*»4  -  */  + 2k^1 


s 


156. 


In  this  case,  the  transferred  power  is  given  by 


14.  RESULTS 


aBa)  v  v  3k 
T7A5*  2 


157. 


Results  have  been  obtained  for  a  steel  plate  of  dimension  l.m  by  0.5  m  and  thickness  6.35  mm.  The 
Young’s  Modulus  of  steel  is  E=  200  GN/m2  and  the  Poisson’s  ratio  is  0.3.  Two  structural  loss  factors 
are  considered  in  the  analysis  10~2  and  10-5.  The  plate  is  excited  by  an  incident  acoustic  waves  with 
pressure  amplitude  10~*  Pa  and  with  angles  of  incidence  of  6  =  45°  and  <p  -  45°  .  Results  are 
generated  with  and  without  the  influence  of  the  fluid  loading.  In  generating  these  results,  both  the 
Kernel  approximation  method  and  the  projection  method  are  used. 


The  results  for  the  inpur  and  transferred  power  are  normalized  with  respect  to  the  incident  power 
on  the  plate,  which  is  defined  by 


nat> 

2g,c 


The  ratio  of  the  input  power,  which  is  the  power  transmitted  to  the  plate,  to  the  incident  power  on 
the  plate  is  a  measure  of  the  transmission  coefficient  T  of  the  plate  (T  =  II  JIT)  and  the  ratio  of  the 
scattered  power  to  the  incident  power  on  the  plate  is  the  portion  of  the  power  scattered  to  the  far  field 
and  is  defined  as  the  reflection  coefficient  R,  (R  =  /7,/P/). 

The  results  for  different  mode  numbers  and  different  levels  of  the  internal  loss  factor  are  shown  in 
figures  (2)  to  (19). 

14.1.  AIR-LOADING 

Figures  (2)  to  (9)  show  the  results  for  the  transmission  and  reflection  coefficients  for  different  modes 
and  loss  factors  using  the  two  alternate  solution  methods  presented  in  sections  (5)  and  (6).  It  can 
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be  observed  from  these  figures  that  in  the  case  of  light  fluid  loading  (air),  the  Kernel  approximation 
method  and  the  projection  method  give  similar  results,  especially  for  high  loss  factors.  At  high 
frequencies  the  two  methods  give  exacdy  the  same  results.  At  low  frequencies  there  are  some 
differences,  especially  in  the  results  for  the  lower  loss  factor  and  these  are  expected  since  in  the 
Kernel  approximation  method  a  short  wavelength  assumption  is  used.  At  low  frequencies  the 
infinite  plate  approximation  does  not  hold.  The  kernel  approximation  method  has  a  tendency  to 
overestimate  the  input  power  to  the  plate.  However  for  the  scattered  power  from  the  plate,  the  results 
given  by  the  kernel  approximation  and  the  projection  method  and  results  obtained  with  the  fluid 
loading  term  neglected  in  the  solution  for  the  response  of  the  plate,  are  the  same.  The  main  reason 
for  this  is  that  the  plate  response  is  hardly  influenced  by  the  light  fluid  loading 

From  these  results  one  can  conclude  that  even  light  fluid-loading  has  some  influence  on  the 
transmission  coefficient  of  the  plate,  at  low  mode  numbers.  The  fluid  loading  increases  the  input 
power  to  the  plate  since  some  power  is  radiated  into  the  acoustic  medium.  The  influence  of  the  fluid 
loading  is  also  dependent  on  the  structural  loss  factor.  From  figures  (6)  and  (7)  it  can  be  observed 
that  away  from  resonance  frequencies,  the  scattered  power  is  independent  of  the  structural  damping. 
On  the  other  hand,  the  input  power  appears  to  increase  with  increased  structural  damping. 

To  show  the  importance  of  acoustic  scattering,  figures  (8)  and  (9)  show  both  the  transmission  and 
the  reflection  coefficients  on  the  same  plot.  The  input  power  into  the  plate  is  always  greater  than 
the  scattered  power  from  the  plate.  This  is  consistent  with  the  requirement  for  the  conservation  of 
energy.  The  input  power  to  the  plate  is  balanced  by  the  scattered  power  and  the  structural  dissipation 
which  is  always  positive,  equation  (86).  At  almost  negligible  loss  factor,  the  input  and  the  scattered 
power  are  almost  equal. 

When  the  structural  damping  is  very  low,  the  resonance  peaks  are  very  sharp  and  no  modal 
interactions  occur  as  the  energy  of  the  system  is  contained  in  very  narrow,  separate  bands  of 
frequencies.  Thus,  the  acoustic  field  at  any  frequency  is  generated  solely  by  the  mode  that  is  resonant 
at  or  near  that  frequency.  The  radiation  damping  is  mainly  due  to  the  acoustic  field  generated  by 
that  mode.  The  total  power  flow  can  thus  be  obtained  from  the  modal  power  flow  for  the  plate. 
When  the  structural  damping  is  large,  the  width  of  the  resonance  peaks  increase,  and  each  mode  is 
excited  over  a  wider  frequency  band.  In  this  case  to  obtain  the  total  power  flow,  a  summation  over 
excited  modes  is  necessary. 

14.2.  WATER-LOADING 

In  the  case  of  heavy  fluid  loading  (water),  figures  (10)  to  (14),  show  that  the  results  obtained  by  the 
Kernel  approximation  method  become  inconsistent  with  those  given  by  the  projection  method.  The 
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kernel  approximation  method  has  a  tendency  to  overestimate  the  input  power  into  the  plate  and 
underestimate  the  scattered  power  from  the  plate. 

Both  methods  predict  a  strong  fluid  loading  effect  on  both  the  transmission  and  reflection 
coefficients.  The  fluid  loading  significantly  modifies  the  plate  mode  shapes.  There  is  a  decrease 
in  the  resonance  frequencies  from  the  in-vacuo  resonances  frequencies.  This  is  because  of  the  added 
mass  entrained  by  the  fluid  and  which  is  significant  in  the  case  of  water.  As  for  the  light  fluid  loading, 
the  structural  loss  factor  has  a  significant  influence  on  the  power  flow  for  the  plate,  figures  (12) 
shows  the  influence  of  the  structural  damping  on  the  transmission  and  reflection  coefficients.  As 
in  the  case  of  light  fluid  loading  (air),  the  scattered  power  from  the  plate  is  always  less  than  the  input 
power  (figure  13),  however  in  the  case  of  the  heavy  fluid  loading  (water),  the  plate  scatters  a 
significant  portion  of  the  power  input  to  the  plate. 

From  all  the  results  for  the  heavy  fluid-loaded  plate,  it  can  be  observed  that  each  resonance 
frequency  is  preceded  by  a  minimum.  The  response  of  the  plate  is  composed  of  an  almost  uniform 
driven  response  as  the  plate  were  infinite  and  an  oscillating  term  due  to  the  edge  effects.  Below  a 
resonance  frequency,  there  is  a  phase  shift  of  it  between  the  driven  response  and  the  oscillating  term 
and  cancellation  occurs,  which  is  responsible  for  the  minimum  value  right  before  the  resonance 
condition.  At  and  above  a  resonance  frequency,  the  two  terms  are  in  phase  and  their  effects  adds. 

At  the  resonance  frequencies,  both  the  transmission  dn  reflection  ratios  have  a  huge  peak.  Near  the 
resonance  frequencies,  the  plate  absorbs  more  power  than  is  directly  incident  upon  it  The  plate 
absorbs  power  from  the  acoustic  waves  incident  on  the  adjacent  baffle.  This  energy  is  then  scattered 
by  the  plate.  At  the  resonance  frequencies,  the  induced  velocity  field  on  the  plate  vibration  becomes 
large  and  increases  the  level  of  scattering.  These  results  are  similar  to  those  found  by  Photiadis  [5] 
for  a  semi-infinite  plate.  Photiadis  suggests  that  the  plate  grabs  energy  incident  on  the  adjacent 
baffle  to  within  a  distance  of  A0/2,  where  A0  is  the  acoustic  wavelength.  In  the  case  considered  here, 
A0  is  larger  than  the  plate  dimensions  at  all  the  considered  frequencies.  Thus  the  energy  absorbed 
by  the  plate  from  the  adjacent  field  is  significant. 

14.3.  COMBINATION  WATER  AND  AIR  LOADING 

Figure  (15)  sho^ ,  i!  ?  transmission  and  reflection  coefficients  computed  using  the  Projection 
Method  when  the  fluid  loading  on  the  excited  side  of  the  plate  is  water  and  the  other  side  the  fluid 
loading  is  air.  The  results  for  the  mixed  fluid  loading  are  compared  to  other  fluid  loading  conditions 
in  figures  (16)  and  (17).  The  scattered  power  is  highest  when  the  plate  is  water  loaded  on  one  side 
and  air  loaded  on  the  other  side.  This  increase  in  the  scattered  power  is  attributed  to  the  way  the  plate 
interacts  with  the  fluid  loading  and  the  excitation. 
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15.  CONCLUSION 


A  solution  has  been  obtained  for  the  dynamic  response  and  acoustic  scattering  from  a 
fluid-loaded  rectangular  plate  with  arbitrary  boundary  conditions  on  two  of  the  parallel  edges  of  the 
plate.  The  solution  takes  into  account  the  finiteness  of  the  problem,  especially  when  dealing  with 
the  scattered  pressure.  The  influence  of  the  finiteness  of  the  scattered  pressure  is  negligible  for  light 
fluid  loading,  however  for  heavy  fluid  loading  it  has  a  significant  influence  on  the  scattered  pressure. 
The  solution  that  is  presented  here  can  be  used  to  analyze  the  response  of  coupled  plate  systems  taken 
into  account  both  the  mechanical  coupling  and  the  acoustic  coupling  through  the  common  acoustic 
medium. 
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Figure  2.  Transmission  coefficient  for  a  simply  supported  plate  with  air  loading,  loss  factor  0.01.  for  different 

modes  and  different  methods  of  solution. - :  Projection  method; - :  Kernel  approximation  method; 

- — No  fluid  loading,  (a)  mode  1;  (b)  mode  2;  (c)  mode  3;  ( d)  mode  4. 


Figure  3.  Reflection  coefficient  for  a  plate  simply  supported  with  air  loading,  loss  factor  0.01,  for  different  modes 

and  different  methods  of  solution. - :  Projection  method; - :  Kernel  approximation  method; 

— - :  No  fluid  loading,  (a)  mode  l;(b)  mode  2;  (c)  mode  3;  (d)  mode  4. 
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Figure  4.  Transmission  coefficient  for  a  simply  supported  plate  with  air  loading,  loss  factor  10~5.for  different 

modes  and  different  methods  of  solution. - :  Projection  method; - :  Kernel  approximation  method; 

— — :  No  fluid  loading,  (a)  mode  1 ;  (b)  mode  2;  (c)  mode  3;  (d)  mode  4. 


Figure  5.  Reflection  coefficient  for  a  simply  supported  plate  with  air  loading,  loss  factor  10~5,for  different  modes 

and  different  methods  of  solution. - :  Projection  method; - :  Kernel  approximation  method; 

- :  No  fluid  loading,  (a)  mode  1;  (b)  mode  2;  (c)  mode  3;  (d)  mode  4. 
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Figure  6.  Transmission  coefficient  computed  by  the  Projection  method  for  a  plate  simply  supported  on  all  edges 

with  air  loading,  for  different  modes  and  different  structural  loss  factors. - :  Loss  factor  of  10~2; - 

Loss  factor  of  10 ~s.  (a)  mode  1;  (b)  mode  2;  (c)  mode  3;  (d)  mode  4. 
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Figure  7.  Reflection  coefficient  computed  by  the  Projection  method  for  a  simply  supported  plate,  for  different 

modes  and  different  structural  loss  factors. - :  Loss  factor  of  10 - :  Loss  factor  of  10 ~5.  (a)  mode 

1;  (b)  mode  2;  (c)  mode  3;  (d)  mode  4. 
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Figure  8.  Transmission  and  Reflection  coefficients  computed  by  the  Projection  method for  a  simply  supported 

plate  with  air  loading,  loss  factor  0.01  .for  different  modes.  ■  Transmission  coefficient; - ; 

Reflection  coefficient,  (a)  mode  1;  (b)  mode(2);  (c)  mode  3  and  (d)  mode  4. 
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Figure  9.  Transmission  and  Reflection  coefficients  computed  by  the  Projection  method  for  a  plate  simply 

supported  with  air  loading,  loss  factor  10~s,for  different  modes. - :  Transmission  coefficient; - : 

Reflection  coefficient,  (a)  mode  1;  (b)  mode(2);  (c)  mode  3  and  (d)  mode  4. 
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Figure  10.  Transmission  coefficient  for  a  simply  supported  plate  with  water  loading,  loss  factor  1  Or2,  for  different 

modes  and  different  methods  of  solution. - :  Projection  method; - :  Kernel  approximation  method; 

- :  No  fluid  loading,  (a)  mode  1;  ( b )  mode  2;  fc)  mode  3;  (d)  mode  4. 


Figure  11.  Reflection  coefficient  for  a  simply  supported  plate  with  water  loading,  loss  factor  lO-^.for  different 

modes  and  different  methods  of  solution.  — — — :  Projection  method; - :  Kernel  approximation  method; 

— — — :  No  fluid  loading,  (a)  mode  1;  (b)  mode  2;  (c)  mode  3;  (d)  mode  4. 
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Figure  12.  Transmission  and  reflection  coefficients  of  mode  1  computed  by  the  Projection  method  for  a  simply 

supported  plate  with  water  loading  for  different  structural  loss  factors.  -  Loss  factor  of  10~2 ; - — ; 

Loss  factor  of  10~ * .  (a)  Transmission  coefficient;  (b)  Reflection  coefficient. 


Figure  13.  Transmission  and  Reflection  coefficients  computed  by  the  Projection  method  for  simply  supported 

plate  with  water  loading,  loss  factor  0.01,  for  different  modes.  —  Transmission  coefficient; - : 

Reflection  coefficient,  (a)  mode  1;  (b)  mode(2);  (c)  mode  3  and  (d)  mode  4. 
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Figure  14.  Transmission  and  Reflection  coefficients  of  mode  1  computed  by  the  Projection  method  for  a  simply 
supported  plate  with  water  loading,  loss  factor  lO"5.  :  Transmission  coefficient;————:  Reflection 

coefficient. 


Figure  15.  Transmission  and  Reflection  coefficients  for  mode  1  for  a  simply  supported  plate  with  waterlair 
loading,  plate  loss  factor  lOr2-. - Transmission  coefficient;— — Reflection  coefficient. 
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Figure  16.  Transmission  coefficient  for  mode  I  for  a  simply  supported  plate  for  different fluid  loadings,  loss  factor 
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Figure  1 7.  Reflection  coefficient for  mode  1  for  a  simply  supported  plate  for  different  fluid  loadings,  loss  factor 
10~2-.  — '  - :  Air  loading;— - :  Water  loading; - :  Water! Air  loading 
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Appendix  A 


rb  ra  '  rb 

P„Jfi)  =  P„Jy)eW»dy  =  §  PSlU.y)emdy  sin (2f£)<fa  A.1 

Jq  Jo  Jo 


Using  the  inverse  Fourier  Transform 


Oao 

. 


PsAa'^’^-VydaW  A.2. 


00  J  —  00 


r  <»  /■  » 


aj*> 


C«8  -  fi')P1:(a',fi')Sm(- a'yki’dP'  A.3. 


00  V  —  00 


where  G(fi)  =  f  e$ydy  =  ~~J6 — "  and  Sm(a)  =  f  sini^—j^dx. 

J  O  Jo 


n  =  » 


H'fcy)  =  2  K(y)sin(-^> 


ft  =  1 


and 


n=  oo 


W(a',F)  =  ^  W„(fi')Sn(a') 


n  =  1 


Substituting  into  equation  (24) 


5  on  _ 

"r,(«  >p  )  =  ;Piw  2^ - y - 


n  =  1 


and 


=  VPPi 

where  n  =  aj2/(2jt2a)  and 


n=  ao  r  * 

xf 

n=  1  J  - 


A.4. 


A.5. 


A.6. 


Wn(fi')G(fi  -  py^ww  A.7. 
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Jyjfi') 


Snia'tfU- 

y' 


a’) 


da' 


A.8. 


Similarly 


wjfi'ycxfi  -  fi'VimJPW 


A.9. 


The  negative  sign  is  introduced  in  equation  (A.9)  because  the  scattered  pressure  waves  are 
propagating  in  the  negative  z  direction. 

In  equations  (A.7)  and  (A.9),  the  terms  of  the  summation  for  which  p  *  m  represent  the  intramodal 
interaction.  This  coupling  has  been  shown  to  be  negligible  when  the  fluid  loading  is  light  and  the 
modal  density  of  the  structure  is  low  [4].  If  the  intramodal  coupling  is  neglected. 


=  2 k2„ 


+  00 

I 


[l  -  (-  l)wcos(aa)j 

(a2  ~  *m)V 


da 


A.  10. 


where  y  =  ±  Jk2  -  (}2  -  a2.  The  sign  of  the  square  root  term  is  chosen  according  to  the  radiation 
condition,  that  is  the  scattered  pressure  is  equal  to  zero  at  infinity. 

If  k  <  ,  y  is  imaginary,  y  —  —  j  y If}2  —  k2  +  a2 

If  k  >  (i,  y  =  [Jk2  —  f}2  —  a2  for  a  s  k2  -  j 32  and  y  -  -  j  Ja2  -  k2  +  fl2  for 
a  "3:  k2  ~  pi2. 

If  k  =  P,  y  is  imaginary,  y  =  -  ja. 

In  equation  (A.  10),  the  integrand  is  an  even  function  of  a  and  therefore 


+  • 


JlJfi)  =  4  kl 


I 


[l  -  (-  l)mcos(aa)] 

(a2  -  4)2y 


da 


A. ll. 


Also  JijP)  is  an  even  function  of  P,  and  therefore  the  integral  in  equation  (A.  10)  can  be  limited 
to  P  s  0.  The  solution  for  different  ranges  of  k  and  P  is  as  follows: 
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Case  L  <.{L 


+  oo 

w-*&\  wh*1 

o 

CassJL=Ai 

+  00 

Jlm(k)  =  4*2/  |  ^rda  A.  13. 

o 

At  zero,  the  integrand  goes  to  infinity  when  m  is  odd  and  equals  to  zero  when  m  is  even.  However 
the  integrand  is  also  an  odd  function  for  odd  values  of  m  and  therefore  the  value  of  the  integral  always 
goes  to  zero  for  odd  values  of  m. 

Caa&-fl  „< 


Jijfi)  -  4 ki 


[l  -  (—  l)“cos(aa)] 


+  oo 


J  (a2- 


(a2  - 


(  1  -  (-  l)mcos(aa) 

da  +  j  I  - - - - Ida 

,  (a*  -  klf 


A.  14. 


Jljfi)  =  4*i 


I 

I 


+  00 


f(VcosX)dX  +  j  f(VchX)dX 


I 


A.15. 


In  equations  (A.  12),  (A.  13)  and  (A.  14),  V  =  J\k2  -  /32|,  and 


[l  -  (-  l)wcos(aa)] 

- - -  A.  16. 

(«2  -  klf 

The  function  f  is  a  smooth  function  and  is  not  singular  at  km.  Fora  =  km)  f{km)  -  a2/(8  m2n2). 

In  equation  (A.  15),  the  substitutions,  x  =  a/V,  and*  =  cos(X)  for  the  first  integral  and  x  =  ch(X) 
for  the  second  integral  have  been  used. 

The  integrands  in  these  integrals  are  smooth  functions  and  therefore  can  be  directly  evaluated 
numerically. 
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Appendix  B 

From  equation  (28),  Axifi)  =  A2(fi)  A3(fl)  and  A4(/3)  =  1 


i=4 

X  =  Atf$AtfVCu  +  Xu]  +  AtfyXu  +  X„  B.1 

i  =  l 

From  equation  (20): 

A3Wu  +  *24  =  -  2 

{sin(AT2ft)[(e -  e#b)K(,Wm)(Kl )  -  («»*  -  ^6)*:(VVm)(  -  AT,)]  - 
sin(A-,i,)[(e -  e>f>b)K(Wm)(t:2)  -  (e>K*  -  ^b)K(Wm)(  -  AT2)]j  „.2 

A30 SVC3A  +  *4A  =  -  ^-{sin(/f2f>)(A:|  +  2*i.) 

[(e  -yjr.*  -  ^b)K(W„XK,)  -  (e**  -  e»b]K(W„X~  X,)]  - 

sin («•,«(*;  +  2^,)[(e-»  -  e^b)K(.W„)(K2)  - 

{e‘K’b  _  ^6)X(W„;(  -  AT2)]}  B.3 

i- 4 

X^A,^)  =  i(sm(^)032  -  K2) 

/  =  1 

[(«'**  _  ^^(I^AT,)  -  (e*»  -  e^W^C-  *,)]  - 
sintA:,*)^2  -  AT2)[(e  -•'**"  -  ^^(IV^X^)  -  (e*-6  -  e*»)KWJ(-  K2)]l  B.4 

Substituting  for  function  K 

i =4  f+ao 

XXiAW  =  -  I  [A,/,.)?')  +  B.5 

i  =  1  J  —  oo 

where 
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K(fi,P')  =  i(</32  -  K]) sm(K2b) 

[C(K,  -  P')(e-iK‘b  -  <**>»)  -  G(-  Kt  -  P')(e‘K‘b  -  e*’4)] 

(fi2  -  AT2) sin(AT,6)f 0(^2  -  P'ie-**  -  e*4)  -  G(-  K2  -  /S')^  -  c^4)] 

B.6 


69 


Appendix  C 

Functions  <pn  are  defined  by  =  d?G(fi)ldnfi,  where  Gift)  =  e$b  -  1/0)3). 

The  recurrence  relation  for  <j)n  may  be  obtained  by  operating  on  the  series,  (3.27). 

00 

G(z  +  0)  =  2  C.1 

n  =  0 


Using  the  definition  for  G(fi) 

eM+V  -  i  =  j(fi  +  j)  y  C.2 

W§  • 

n-0 


Differentiating  this  equation  with  respect  to  z  and  then  with  respect  to  (i 

=  £  j£[(*  +  2)^„+103)  +  p<pn+2(fi)}  C.3 

n  =  0 

where  use  was  made  of  the  relationship,  <pm+  j(|8)  =  chpm(fi)/ (dfi) 

Substituting  for  the  exponential  term 

<pn-m  =  ^((n  + 1 ~ +  C4 

Given  <p0  and  <j>x  any  other  <pn  can  be  computed.  For  high  value  of  n  and  small  values  of  /?,  <pn(fi) 
is  very  small  and  numerical  problems  occur.  A  better  way  to  compute  these  functions  is  to  use 
Miller’s  algorithm  applying  the  recurrence  downward  [4  page  142  and  28].  The  method  of 
computing  <pn(fi)  consists  of  starting  from  some  arbitrary  value  <Pm0)  —  e  with  M  »  n  and  e  a 
small  number,  and  making  use  of  the  recurrence  downward.  All  the  <pn(fi)  thus  computed  have  a 
common  normalizing  factor  given  by 

4>n(°)  =  JTjW  C.5 

b  =  Yhr-t'W  C6 

ft  • 

0 
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In  computing  <  <pn,  <f>t  > , 


Appendix  D 


<  <Pn*<Pi  >  - 


I 


<PnW,  *  W 


where  <p„(fi)  =  dnG(j$)/dfin.  The  function  G  can  be  expressed  as  an  integral 

b 


and  therefore 


<Pntf)  = 


D.l 


D.3 


b  b  +« 

<Pn,<Pi  >  =  j  J(/yi),,(/y2)/  | 

0  0  L-® 


.n  +  l 


dyidy2  =  2 Jib(- 


D.4 


since 


+  «0 


I 


—  00 


eiP^-y^dp  =  2nd(yl  -  y2 ) 


D.5 


In  computing  <  rp„,<f>i  > 


<  V>n,<Pi  >  = 


f 


V-»<fl0/  *  (flop 


D.6 


and  substituting  for  front  equation  (61)  and  interchanging  the  order  of  integration. 
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+  00 

M,  >  =  J  [A,y,_03->  + 

—  00 

f  +  * 

?,03'J  -  j  lG(fi  -  F)  +  K<P,P')W  */  (fiw 


A  contour  integral  can  be  used  to  evaluate  eq’  3  ion  (D.8).  Substituting  for  <pt  *  (fi)  from  equation 
(D.3)  and  substituting  for  G(fi  -  0')  from  equation  (D.2). 

0,0')  =  (-  jj  J>'  [J  ^jG03  -  P')e-fod$  +  j  ^K00')e-^dfi  dy 


J  —  oo 


yf  +* 


*r  + 00 


e^Cy.-y) 

A0) 


<#kyi  + 


0  L-» 


e^Cyi-y) 

A09) 


>  L- 


yr  +  * 

f  f  e^(y»-y)^a  .  _  «r  e^»y  -  g-ff'y  ,  -  e~^>  1 

J  J  A(jS)  ^  [(AT,  +  /9'VW,)  (*2  +  /3>'(^2)J 

0  00 


61"  +  00 

f  f  eJtoi-y)  JO  .  „  r eM'-We-JK#  ■ 

I  J  ».  -  (if,  _  ^ 


-jK,y  -  e-p,  ej(K2-py>e-jK#  -  e-j/l'y 
P')A'(Kl)  (K2  ~  P')A'(K2) 


72 


where 


Htf')  =  -  2kp  sin(AT26)[G(  —  AT,  -  ?)  -  G(A",  -  /S')] 

W2(/S')  =  -  2itpSm(A"1i)[G(-  AT2  -  £')  -  G(K2  -  j9’>]  D.14 

//3</S')  =  -  2/tJ  sin(it2i>)[e  ~iK'bG(K j  -  /S')  -  e^'*G(-  -  /S')] 

Ht0‘)  =  -  sin(AT (6)[e  ~ik^G(K2  -  /S')  -  eiK’bG(  -  K2  -  /3')] 

<pfK),  for  K~  ±  Kx ,  and  ±  AT2 ,  can  be  evaluated  using  equation  (D.3). 

In  computing  <  WQm,<pt  > 


+  ao 


<W^/>=  j 
—  00 

Using  equation  (D.3) 

wajfiyp  *, 

b 

r 

r  +  00 

r 

<  »'o..0/>=  <-/y)' 

Wo_0)e  "*'43 

0 

J  ! 

-00 

Inverse  Fourier  transforming  the  integral  between  the  brackets 

b 

r 

<  W0m,<pt  >  =  2 ji 

(-  jy)lwmSy)dy 

0 
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Structural  Intensity  Measurements 

There  are  a  variety  of  wave  modes  which  can  propagate  in  a  vibrating  structure.  Some  waves 
can  propagate  almost  unaffected  past  energy  dissipation  or  isolation  mechanisms  meant  for  other 
wave  modes.  This  has  the  effect  of  ‘short  circuiting’  the  reduction  mechanism(s)  for  vibrational 
energy  propagation.  One  example  is  the  case  of  in-plane  and  out-of-plane  waves  propagating  in 
a  structure.  The  in-plane  waves  propagate  past  the  isolators  and  damping  layers  used  to  reduce  the 
propagation  of  bending  waves.  Furthermore,  some  wave  types  are  more  efficient  acoustic  radiators 
than  others,  and  can  be  the  dominant  sources  of  radiation  from  a  vibrating  structure.  For  example, 
subsonically  propagating  out-of-plane  waves  radiate  inefficiently  in  water,  in  contrast  with  the 
supersonic  in-plane  waves  which  are  efficient  radiators. 

The  scattering  of  propagating  waves  by  discontinuities  in  their  paths  can  alter  their  relative 
energy  content.  For  instance,  the  energy  of  incident  bending  waves  can  be  redistributed  into  in-plane 
longitudinal  and/or  in-plane  shear  waves  when  the  waves  encounter  discontinuities.  The  scattering 
processes)  can  thus  lead  to  an  increase  in  the  acoustic  radiation  from  a  structure.  That  is,  the  radiated 
noise  is  a  function  of  the(se)  scattering  process(es).  Thus,  the  relationship  between  the  nature  of  the 
discontinuities  and  the  various  wave  scattering  mechanisms  needs  to  be  better  understood.  This 
would  enable  the  governing  parameters  of  the  power  exchange  between  various  wave  modes  to  be 
quantitatively  identified. 

A  useful  tool  to  describe  the  dynamic  behavior  of  structures,  is  the  concept  of  structural 
intensity.  In  general,  the  instantaneous  structural  intensity  (or  Poynting)  vector  at  a  point  with 
position  vector?,  in  an  isotropic,  lossless,  elastic  medium,  is  defined  [1]  as, 

J(r,t)=  -W!!dl.o(r,t)  17. 

where  J (P,  t)  is  the  displacement  vector  and  a  (7,  t)  is  the  stress  tensor 

The  magnitude  of  the  structural  intensity  vector  is  thus  the  power  associated  with  the  waves 
propagating  in  the  elastic  medium,  and  the  direction  of  the  vector  is  the  direction  of  propagation  of 
this  power.  For  the  more  realistic  case  of  losses  existing  in  the  elastic  medium,  there  will  be  a  phase 
lag  ‘<J>’  between  the  stress  and  displacement.  Thus,  material  damping  is  accounted  for  by 
incorporating  this  phase  lag  in  the  time  dependent  parts  of  the  expression  for  the  structural  intensity. 

It  is  of  interest  to  consider  two  simplifications  of  this  definition  [2].  Firstly,  the  instantaneous 
intensity  can  be  time-integrated  (for  transient  waves)  or  time-averaged  (for  steady-state  waves). 
This  yields  a  frequency-domain  (using  Parseval’s  thee  ri  description  of  the  net  energy  flow  or 
time-averaged  power  flow  respectively,  7  (r,  to),  whic,  re  useful  for  source/sink  localization. 
The  second  simplification  involves  specialization  to  two-dimensional  structures  such  as  plates  and 
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shells.  In  this  case,  it  is  usual  to  assume  that  the  displacement  vector  s  (r,  t)  can  be  expanded  in  a 
series  in  the  thickness  direction,  about  the  middle  surface  of  the  cross-section.  Substitution  in  the 
expression  for  the  intensity  vector,  and  integration  over  the  thickness  of  the  structure,  leads  to  the 
simplification  of  the  intensity  as  an  intensity  resultant  per  unit  length. 


It  can  also  be  shown  that  upon  rearrangement  of  the  expressions,  one  obtains  the  intensity 
resultant  vector  as  a  product  of  the  stress  resultant  tensor  and  the  corresponding  velocity  vector.  It 
is  thus  possible  to  decompose  the  total  structural  intensity  resultant  vector  7s,  into  the  (scalar  )  power 
flow  contributions  associated  with  each  wave-type  and/or  direction.  For  two-dimensional 
structures,  these  contributions  would  come  from  the  in-plane  shear,  in-plane  longitudinal  and 
out-of-plane  waves.  Each  of  these  wave-types  has  different  dispersion  characteristics  with 


corresponding  propagation  wavenumbers,  for  the  in-plane  shear  waves,  ks  =  Inf/UG/g j,  for  the 
in-plane  longitudinal  waves,  kL  =  2nf/{jE/o}  and  for  the  out-of-plane  (transverse)  waves 


J2rf  (Eh2/(e)) 


-1/4 


The  component  wavenumbers  each  have  a  specific  dependence  on  frequency  and  the 
medium’s  inertial  &  elastic  properties.  Thus,  a  wavenumber  decomposition  of  the  structural 
intensity  would  enable  identification  of  the  contributions  from  the  various  propagating  waves.  This 
decomposition  can  lead  to  a  better  understanding  of  the  mechanisms  of  scattering  by  discontinuities, 
modal  energy  conversion  &  distribution  and  modal  coupling.  This  understanding  is  of  importance 
in  the  investigation  of  the  dynamic  response  of,  and  acoustic  radiation  from,  fluid-loaded  structures 
consisting  of  coupled  compk  ibstructures. 

This  paper  deals  win  theoretical  analysis  of  structural  intensity  measurement  schemes 
for  thick  plate  structures.  Firstly,  the  frequency-wavenumber  domain  structural  intensity  is  defined. 
Then,  the  various  intensity  components  for  the  case  of  thick  structures  are  described.  Next,  a  scheme 
is  described  for  measuring  the  frequency-wavenumber  domain  structural  intensity,  which  uses  finite 
difference  approximations  for  the  various  spatial  derivatives  required.  The  proposed  scheme’s 
performance  is  considered  for  the  practical  case  of  relatively  small  data  sets,  and  for  various  levels 
of  Signal-to-Noise  Ratio  (SNR).  To  study  the  limits  of  applicability  of  the  structural  intensity 
measurement  scheme,  it’s  performance  is  characterized.  Finally,  the  results  of  experiments  using  the 
proposed  measurement  scheme  are  also  presented  in  order  to  illustrate  the  various  issues  involved. 


Frequency  -  Wavenumber  Analysis 

From  the  definition  of  structural  intensity  in  equation  (1),  the  intensity  resultant  over  a 
surface  of  area  5  is  defined  at  an  instant  of  time  t  as, 
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18. 


r)  =  ^  I  F(7,  t).v(r,  t).nr.dS 
Js 

where  T  is  the  position  vector,  and  nr  is  the  unit  normal  to  S  and  T,  and  F(r,  t)  and  V(r,  t)  are  the 
generalized  force  and  velocity  vectors. 

The  time-averaged  (time -integrated)  structural  intensity  for  the  case  of  real,  steady-state 
(transient)  time-series,  is  given  from  an  application  of  the  power  theorem  [1]  by, 

F(7,a>)  -  Real  [ F(r,o)).V  (7,  to)]  19. 

where  F(f,  to)  and  V(r,a>)  are  respectively  the  force  and  velocity  vectors  in  the  transformed 
frequency  -  space  domain  and  *  represents  a  complex  conjugate. 

For  the  analogous  case  of  spatial  transformation  to  the  wavenumber  domain,  it  is  useful  to 
consider  a  general  vector  spatial  field,  where  the  position  of  a  point  is  defined  by  a  vector  r.  The 
Held  consists  of  a  generalized  force  F  and  a  generalized  velocity  v  which  vary  with  both  space  and 
time.  The  cross-correlation  Rpv  of  the  force  and  the  velocity  is  defined  using  the  expectation 
operator  E{  }  as, 

R^  (7,1 1,  X)  =  E  {  F(7,  t).v(7  +  l  t  +  X)  }  20. 

The  force-velocity  cross-spectrum  Sp  v  ,  given  by  the  Fourier  transformation  of  the 
cross-correlation  over  time  t  and  space  r  ,  gives  the  frequency-wavenumber  domain  structural 
intensity.  For  a  stationary  field,  the  cross-correlation  Rp  v  (and  every  higher  statistical  moment  in 
time)  is  independent  of  the  absolute  observation  time  t ,  and  depends  only  on  the  observation  time 
difference  T.  Similarly,  for  a  homogeneous  field,  the  cross-correlation  (and  every  higher  statistical 
moment  in  space)  is  independent  of  the  absolute  spatial  position  r ,  and  depends  only  on  the  spatial 
separation  vector^.  As  steady-state  phenomena  are  of  greater  interest  in  the  present  study,  temporal 
stadonarity  is  assumed.  The  structural  discontinuities  are  of  interest  however,  and  consequently, 
spatial  homogeneity  is  not  assumed.  For  the  homogeneous  case,  the  cross  spectrum  has  a 
dependency  only  on  the  wavenumber  H  and  the  frequency  co,  thus  yielding  Spv(k  (0).  However, 
for  the  inhomogeneous  case,  there  is  an  additional  dependency  on  the  absolute  spatial  position  r , 
apart  from  the  wavenumber  k  and  the  frequency  (ft  which  yields  Spv(kt  (0,  r).  This  additional 
dependency  on  the  absolute  position  is  resolved  by  spatial  averaging  the  cross-spectrum  over 
sections  of  the  domain  of  interest  for  which  the  spatial  field  can  be  considered  homogeneous.  This 
gives  the  frequency-wavenumber  domain  cross-spectrum.  Thus,  taking  a  spatial  average  (integral) 


over  die  domain  of  the  structure,  the  structural  intensity  in  the  frequency-wavenumber  domain  is, 
again  from  the  power  theorem  [1]  given  by, 


P(k,(D )  =  Rea^F(k,a)).v  *  (k,Q))^  21. 

If  a  section  in  the  spatial  domain  cannot  be  identified  for  which  the  response  field  can  be 
considered  homogeneous,  then  the  frequency-wavenumber  spectrum  of  the  structural  intensity 
would  be  given  by  the  convolution  integral, 


F(lw)  =  Real  I  ^F(V,oj).V*  (k  -  V,o))^dV 
J  -  00 


22. 


Structural  Intensity  Components  and  Their  Estimation. 

The  intensity  components  to  be  considered  in  a  structure  of  interest  depend  upon  the 
approximation  to  exact  elasticity  theory  used  to  describe  the  structure.  In  thick  plate  structures  where 
the  out-of-plane  motion  is  significantly  affected  by  the  transverse  shear  and  rotary  inertia,  the 
Timoshenko-Mindlin  theory  [2]  is  appropriate.  Based  on  this  theory,  there  are  waves  associated  with 
the  shear,  bending,  and  twisting  moments  for  the  flexural  motion.  The  corresponding  structural 
intensity  components  associated  with  the  out-of-plane  waves  in  the  Jt-direction  are  listed  below. 

The  out-of-plane  shear  force  component  is 

-  \  Real  [G',(f  +  y-,)  23. 

The  out-of-plane  bending  moment  component  is 


The  out-of-plane  twisting  moment  component  is 


lix  =  \  Real 


(1  -  v)D 


25. 


where  v,  D,  G\  E  and  h  are  respectively  the  poisson’s  ratio,  the  flexural  rigidity,  the  corrected  shear 
modulus,  the  modulus  of  elasticity  and  the  thickness,  u,  v,  and  w  are  the  x,  y,  z  direction 
displacements  and  ipx  =  dw/dx  and  rpy  =  dw/dy. 
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The  expressions  for  the  y-direction  components  are  similar  but  with  x  and  y  interchanged. 
In  structures  where  the  coupling  between  the  in-plane  extensional  and  shearing  modes  significantly 
affects  the  in-plane  motion,  an  approach  along  the  lines  of  Kane  &  Mindlin  [3],  Mindlin  and  Medick 
[41  and  Dubbelday  [5]  needs  to  be  considered.  The  present  analysis  will  consider  the  approach  of 
Kane  &  Mindlin  for  the  in-plane  motion,  in  conjunction  with  the  Mindlin  thick-plate  formulation 
for  out-of-plane  motion.  The  corresponding  x-direction  components  of  structural  intensity, 
associated  with  the  in-plane  waves  are  listed  below 

The  in-plane  longitudinal  force  component  is 


Inx  =  \  Real 
The  in-plane  shear  force  component  is 

Isx  =  \  Real 
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27. 


The  expressions  for  the  y-direction  components  are  similar. 

The  hybrid  approach  used  above  has  the  implication  that  there  is  no  direct  coupling  between 
the  in-plane  and  the  transverse  motion.  However,  there  is  an  indirect  relationship  between  the 
flexural  and  extensional  displacements  [6]  that  comes  from  the  generalized  plane  stress  assumption 
for  the  in-plane  motion.  For  plates  lying  in  the  x,  y  plane,  this  implies  that 


ez=  ~  +  ey) 


28. 


Substituting  the  plate  displacements  for  the  strains  in  equation  (12),  the  out-of-plane  motion 
contribution  from  in-plane  waves  in  each  orthogonal  direction  is 


w 


in  -plane,  i 


vh 

2(1  -  V) 
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where  i  =x  or  y.  Thus,  equations  (7)  to  (11)  can  all  be  obtained  in  terms  of  the  out-of-plane 
displacement  V. 

3.  Estimation  of  the  Spatial  Derivatives 

One  method  of  obtaining  the  spatial  derivative  (and  consequently  intensity  )  in  the 
wavenumber  domain  follows  from  Fourier  transform  properties  [1].  Thus,  the  first  order 
derivative(s)  would  be  given  by 
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Vw(7r,  at)  =  $~l[(j]c).w(ic,a))] 


30. 


where  wife,  at)  is  the  spatial  Fourier  transform  (?)  of  a>) 

The  higher  order  derivatives  are  obtained  by  further  application  of  eq.(12).  In  die  presence 
of  noise,  these  estimates  have  been  shown  by  Pate  et  al  [7]  to  degrade  progressively  as  the  order  of 
the  derivadve(s)  increases,  due  to  the  uncertainty  in  wavenumber  components  ^  anci  ky  Instead  of 
this,  finite  difference  approximations  for  the  derivatives  can  be  used,  as  suggested  by  Pavi$  [8]. 
Thus,  we  have  for  the  the  first  order  spatial  derivative(s)  of  the  displacement  in  the  ‘i ’  direction  using 
transducers  spaced  a  distance  ‘A,’  apart, 

dw(r,(o )  _  w(7  +  (o)  -  w(r  -  Au  to) 

dr,.  *  2Ai 

The  higher  order  derivatives  required  to  compute  the  various  generalized  forces  and 
velocities  can  thus  be  similarly  approximated  using  finite  differences,  in  terms  of  the  signals  from 
an  array  of  transverse  displacement  transducers.  The  highest  order  of  spatial  derivatives  in  equations 
(7)  to  (11)  for  thick  plates  is  2.  Consequently,  the  number  of  measurements  required  in  the  array 
are  9,  to  enable  approximation  of  all  the  partial  and  mixed  spatial  derivatives.  Thus  the  array 
geometry  shown  in  figure  (1),  can  be  used  to  estimate  all  of  the  spatial  derivatives  up  to  order  two, 
thus  enabling  measurement  of  in-plane  wave  intensity  as  well  as  the  thick-plate  out-of-plane  wave 
intensity. 


©  ©  © 
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Figure  1.  Array  geometry  with  9  channels  for  measuring  structural  intensity  using  a  finite 

difference  approximation  for  the  derivatives. 

From  equations  (7)  to  (11),  the  structural  intensity  components  can  be  obtained  by  either 
evaluating  convolution  integrals  of  the  form  of  equation  (6).  The  elements  in  the  convolution 
integral  are  obtained  directiy  from  the  spatial  transforms  of  the  response  and  its  derivatives. 
Alternatively,  if  the  response  field  over  a  section  of  the  structure  can  be  considered  homogeneous, 
as  would  be  the  case  under  resonant  behavior,  then  using  equation  (5)  and  substituting  the  finite 
difference  estimates  for  the  required  partial  derivatives  leads  to  expressions  for  the  frequency 
wavenumber  intensity  components  T  associated  with  the  different  wave-types  and/or  directions  as 
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a  linear  combination  of  frequency-wavenumber  cross-spectra  \S/m’  of  the  signals  from  the  array 
transducer  channels  7’  and  'm\  as  follows , 


Pifco))  =  I  alm.Slm(k, (l))  32. 

lm 

where  'dim'  are  the  coefficients  of  the  cross-spectra  ‘5^,’  introduced  from  the  finite  difference 
approximation  and  the  relationships  between  forces  and  measured  responses.  /  and  m  can  take 
values  from  1  to  9. 

This  procedure  is  graphically  represented  in  figure  (2). 


Figure  2.  Measurement  procedure  for  the  frequency-wavenumber  domain  structural  intensity. 
The  wavenumber  transform  is  taken  over  a  spatial  region  of  the  structure  (length  or  area). 

Whether  using  the  space  averaged  result  or  the  convolution  approach,  the  structural  intensity 
estimates  are  only  subject  to  finite-difference  errors,  provided  phase  and  magnitude  errors  in  the 
measurement  have  been  eliminated  or  minimized.  The  finite  difference  errors  can  be  limited  by 
choosing  the  array  dimensions  appropriately  [8,9].  The  expressions  for  the  ^-components  of  the 
finite-difference  based  approximations  of  the  piece  wise  space  averaged  structural  intensity 
components  associated  with  the  various  wave-types,  in  terms  of  the  velocity  cross-spectra  are  listed 
below. 

The  out-of-plane  shear  force  component  is 

‘v=\  Imag  -  S12)]  33. 

The  out-of-plane  bending  moment  component  is 
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34. 


Im*  38  2  Imag  [4^3  {(^33  ~  sn  -  2^23  +  2^21) 

+  v  (583  —  S8i  +  S53  —  S51  —  2523  +  252i)  }  ] 

The  out-of-plane  twisting  moment  component  is 

I*  -  j  Imag  ^  ^  (598  -  S95  +  5^  -  S45  -  578  +  S75  -  +  S65)j  35. 

The  in-plane  longitudinal  component  is 

'«  =  2  [2dJ(f-  -  S“  *>  +  •£’’«  -  *"«))]  36. 

The  in-plane  shear  component  is 

'«  “  5  [4—^^“^  -  +  5”32  -  «"„)]  37. 

where  S;j  is  the  velocity  cross-spectrum  between  channels  ‘i’  and  ‘j’.  The  superscripts  refer  to  the 
in-plane  velocity  components.  The  y-direction  expressions  are  similar.  Each  of  the  cross  spectra 
are  functions  of  frequency  and  wavenumber. 

Measurement  Limits  for  Structural  Intensity. 

The  transformation  of  the  frequency-space  domain  structural  intensity  to  the 
frequency-wavenumber  domain  involves  the  wavenumber  spectrum  estimation  from  the 
space-domain  data.  There  are  a  number  of  methods  available  for  this,  including  Fast  Fourier 
Transform  (FFT),  Auto  Regressive  (AR),  Moving  Average  (MA),  Auto  Regressive  Moving  Average 
(ARMA),  Pisarenko  Harmonic  Decomposition  (PHD),  Maximum  Entropy  Method  (MEM)  and 
Maximum  Likelihood  Method  (MLM). 

The  usual  method  for  estimating  the  spectrum  of  discretely  sampled  deterministic  and 
stochastic  processes  is  based  on  the  FFT.  This  approach  has  the  advantage  of  being  computationally 
efficient  and  is  valid  for  a  large  class  of  signal  processes.  The  most  prominent  limitation  of  the  FFT 
approach  is  that  of  relatively  low  spectral  resolution.  A  second  limitation  inherent  in  the  method  is 
the  energy  leakage  in  the  spectral  domain  from  the  main  lobe  of  the  spectral  response  to  the 
sidelobes.  These  two  limitations  have  an  especially  detrimental  effect  when  only  short  data  record 
lengths  are  available. 

The  windowing  of  data  is  the  fundamental  factor  that  restricts  the  FFT’s  frequency 
resolution.  This  is  a  consequence  of  the  implicit  assumption  that  the  finite  length  of  measured  data 
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is  obtained  by  multiplying  an  infinite  length  of  the  sampled  data  sequence  by  a  window  function 
which  is  zero  outside  the  measurement  interval.  This  multiplication  of  the  actual  time-series  by  a 
windowing  function  implies  that  the  overall  transform  is  the  convolution  of  the  desired  transform 
of  the  actual  time-series  with  the  transform  of  the  windowing  function.  This  convolution  results  in 
the  spreading  of  a  narrow-bandwidth  signal  into  adjacent  frequencies.  This  is  the  leakage 
phenomenon,  which  distorts  the  spectrum  estimate  from  the  desired  true  spectrum. 

The  most  narrow  spectral  response  of  the  transform  estimate  is  limited  to  that  of  the 
main-lobe  width  of  the  window’s  transform,  and  is  independent  of  the  data.  It  is  proportional  to  the 
reciprocal  of  the  length  of  data  being  transformed.  The  type  of  window  function  chosen  depends  on 
the  desired  compromise  between  reducing  leakage  and  increasing  resolution.  This  trade-off  occurs 
because  when  windows  with  tapered  ends  are  used  to  obtain  lower  sidelobes,  the  main  lobe  width 
increases,  reducing  resolution. 

The  various  alternative  methods  listed  previously  have  been  proposed  to  alleviate  the 
inherent  limitations  of  the  FFT.  These  methods  are  often  more  restricted  in  the  class  of  signals  to 
which  they  are  applicable.  They  can  provide  a  substantial  improvement  in  the  spectral  resolution 
and  signal  detectability.  These  performance  advantages  however,  are  strongly  dependent  on  the 
Signal  to  Noise  Ratio  (SNR).  It  has  been  shown  in  previous  studies  [10]  that  for  low  enough  SNRs, 
the  alternative  spectral  estimates  often  provide  no  improvement  over  conventional  FFT-based 
estimates.  They  have  the  further  disadvantage  of  higher  computational  requirements.  The  structural 
intensity  measurements  can  involve  SNR  values  as  low  as  -10  dB.  At  such  low  values  of  SNR,  it 
is  clearly  advantageous  to  consider  an  FFT  implementation  of  the  space-wavenumber 
transformation.  This  provides  a  conservative  estimate  of  the  measurement  limits,  which  can  then 
be  improved  using  alternative  methods  if  necessary. 

Of  the  two  FFT-based  spectral  estimation  techniques  available,  this  analysis  will  consider 
the  periodogram  method  in  preference  to  the  autocorrelation  method.  This  is  because  the 
autocorrelation  estimate  can  sometimes  yield  negative  values  for  the  PSD,  which  is  physically 
unrealistic.  Thus,  when  measured  data  is  available  only  for  samples  Wo,  Wj,  vvy, . . .  ,w^-i ,  equally 
spaced  Ax  apart,  the  estimate  of  the  spectrum  is  given  by. 


WPER  (*)  =  E 


N-  1 

Ax  2 
n  =  0 


w„  e 


-(j  2jz  knAx) 
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which  is  defined  for  the  spectral  interval  -  1/(2  Ax)  <  k  ^  1/(2  Ax).  It  is  to  be  noted  that  the 
expectation  operator  still  needs  to  be  applied.  The  use  of  the  FFT  thus  permits  evaluation  of  the 


WPER  (it)  ata  discrete  set  of  N  equally  spaced  spectral  points  km  =  m  Ak  where  m  =  0, 1, ...,  N-l 
and  Ak  =  1  /(N  Ax). 

For  the  frequency-wavenumber  domain  structural  intensity  measurement  using  finite  and 
small-length  data  sets,  the  wavenumber  resolution  is  the  limiting  factor.  The  dispersion  curves  for 
the  various  wave-types  approach  each  other  as  the  frequency  approaches  zero.  Consequently,  there 
is  a  lower  wavenumber  limit,  knHn,  below  which  the  wavenumber  separation  is  lower  than  the 
available  wavenumber  resolution  Ak  .  The  frequency  corresponding  to  this  wavenumber  is  the 
lower  frequency  limit  Below  this  lower  frequency  bound,  the  FFT-based  measurement  cannot 
distinguish  between  the  dispersion  curves  associated  with  the  various  wave-types.  The  upper 
frequency  bound  is  the  frequency  corresponding  to  the  maximum  wavenumber  that  can  be 
measured,  and  is  given  by  the  Nyquist  criterion  to  be  —  nj (Ax)-  Above  this  wavenumber, 
aliasing  phenomena  will  cause  errors  in  the  estimation  of  the  spectrum. 

Experimental  Analysis 

Experiments  were  performed  using  the  proposed  approach,  on  a  Plexiglass  model  of  a 
semi-infinite  beam,  0.0254m  thick,  1.17m  long  and  0.05m  wide.  One  end  of  the  beam  was 
embedded  in  sand  to  provide  an  anechoic  termination,  as  shown  in  figure  (3).  The  excitation  was 
provided  using  an  electrodynamic  shaker  attached  to  the  beam  by  direct  threading  at  the  free  end. 


Figure  3.  Experimental  set-up  for  frequency-wavenumber  intensity  measurements 

The  frequency-wavenumber  spectrum  was  obtained  by  measuring  the  response  at  179 
locations  spaced  0.012m  apart  The  sample  record  was  zero-padded  to  256  points,  and 
spatial-transformed.  The  limiting  wavenumber  is  given  by 


k 


max 


<  a  _  n 
A  0.012 


261.7m~l 
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where  A  is  the  measurement  spacing.  The  beam-response  measurement  was  performed  using  a 
Radio  Shack  stereophone  cartridge  as  a  velocity  transducer.  The  model  used  had  an  operational 
upper  frequency  limit  of  18  kHz.  The  excitadon  was  consequendy  limited  to  20  kHz.  The  frequency 
response  estimate  at  each  measurement  locadon  on  the  beam  was  obtained  as  a  cross-spectrum 
between  the  acceleradon  at  that  location  and  the  input  load.  The  spatial  transformed  results  for  the 
acceleration  and  its  spatial  derivatives  (obtained  using  finite  differences)  were  then  used  to  compute 
the  structural  intensity  in  the  frequency-wavenumber  domain.. 

The  approach  used  is  that  based  on  a  homogeneous  spatial  field.  That  is  the  structural 
intensity  components  are  evaluated  using  expressions  of  the  form  of  equations  (17)  to  (21).  This 
assumption  is  not  quite  correct  in  the  case  of  the  semi-infinite  beam,  and  as  will  be  shown  later  the 
results  reflect  this.  The  results  of  the  analysis  are  shown  in  figure  (4).  The  two  intensity  components 
observed  correspond  to  the  in-plane  waves  (the  linear  frequency-wavenumber  relation)  and  the 
out-of-plane  waves  (the  non-linear  frequency-wavenumber  relation).  The  horizontal  axis 
corresponds  to  positive  and  negative  wavenumbers  bounded  by  the  maximum  value  computed 
above.  The  vertical  axis  corresponds  to  the  frequency  axis  with  the  18kHz  upper  limit  mentioned 
previously. 


Figure  4.  Experimentally  measured  frequency-wavenumber-domain  intensity 
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Careful  examination  of  the  results  in  figure  (4)  reveals  that  while  the  general  trends  of  the 
two  components  of  intensity  are  as  expected  followed,  the  curves  show  some  features  which  cannot 
be  attributed  to  physical  phenomena.  For  example,  upon  expansion  of  the  wavenumber  region 
containing  the  two  curves  it  is  observed  that  one  of  the  curves  does  not  start  from  the  origin.  These 
inaccuracies  are  attributed  to  the  fact  that  for  a  semi-infinite  beam,  the  response  and  structural 
intensity  are  monotonically  decreasing  away  from  the  input  Therefore  the  approximation  of  a 
piecewise  homogeneous  spatial  field  is  not  appropriate.  In  this  case  a  convolution  approach  to 
evaluate  the  structural  intensity  components  from  the  estimates  of  the  response  and  its  derivatives 
should  be  used. 

Conclusion 

A  method  for  discriminating  between  the  components  of  the  structural  intensity  propagating 
through  a  structure  has  been  developed.  The  dependency  of  the  measurement  method  and  the 
evaluation  of  the  structural  intensity  as  a  frequency-wavenumber  function  for  both  homogeneous 
and  non-homogeneous  response  fields  have  been  explored.  The  effects  of  expected  noise  and 
damping  levels,  as  well  as  the  relatively  short  data  sets,  on  the  measurements  have  been  considered 
to  obtain  an  understanding  of  the  governing  factors  that  control  the  quality  and  resolution  of  the 
results. 

The  proposed  approach  has  been  implemented  on  a  semi-infinite  beam,  using  the 
assumption  of  a  homogeneous  spatial  response  field.  In  this  case,  while  the  expected  general  trends 
can  be  observed  in  the  results,  there  are  inaccuracies  which  can  only  be  attributed  to  the 
homogeneous  field  assumption.  Analysis  of  the  data  without  using  this  assumption  will  be  explored 
in  the  future  together  with  the;  analysis  of  the  structural  intensity  for  the  beam  if  a  structural 
discontinuity  is  introduced  r  ome  distance  away  from  the  excitation  point 
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